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In this paper we find that in the thermodynamic hmit and for the the ground-state normal- 
ordered ID Hubbard model the wave function of excited states contained in few-electron excitations 
O ' factorizes for all values of the on-site Coulombian repulsion U. This factorization results from 

the non-interacting character of the pseudofermions whose occupancy configurations describe these 
' excited states. Our study includes the introduction of the pseudoparticle - pseudofermion unitary 

transformation and of an operator algebra for both the pseudoparticles and the pseudofermions. 
The pseudofermion description takes into account the relationship between the rotated electrons 
and the holons, spinous, and cO pseudoparticles. [Rotated electrons are related to the electrons 
by a canonical transformation.] As the corresponding pseudoparticles, the cv pseudofermions (and 
sv pseudofermions) are 77-spin zero 2z^-holon composite quantum objects (and spin zero 2z^-spinon 
5— ( . composite quantum objects) where v = 1, 2, .... The pseudofermions are non interacting and thus 

j/j ' have no residual interactions, in contrast to the corresponding pseudoparticles, whose statistics 

we classify according to the generalized Pauli principle. The physics behind the invariance of the 
. pseudoparticles under the above transformations for specific values of the bare momentum is also 

studied and discussed. 

' PACS numbers: 71.10.Pm, 03.65.-w, 71.27.+a, 72.15.Nj 
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O i I. INTRODUCTION 

CNJ , Recently there has been a renewed experimental interest in the exotic one-electron and two-electron spectral prop- 
■ erties of quasi-lD materials 0, S IE IE IS- Some of these experimental studies observed unusual finite- 
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energy /frequency spectral properties, which are far from being well understood. For low values of the energy, the 
microscopic electronic properties of these materials are usually described by systems of coupled chains. For finite values 
of the energy larger than the transfer integrals for electron hopping between the chains^he one-electron (ID) Hubbard 
model is expected to provide a good description of the physics of these materials |E IE IE S & iS] ■ This is confirmed 
ff^ ' by the recent quantitative studies of Refs. P, 0| . Similar unusual spectral properties observed in two-dimensional 
, (2D) high- Tc superconductors could result from effective quasi-lD charge and spin transport HTl IT^ IT?. Ti*. TF] . 
' However, the non-perturbative nature of the ID Hubbard model implies that the electronic creation and annihilation 
^ ] operators do not provide a suitable operational description for the study of the finite-energy spectral properties. Thus, 
G ■ the first step for the study of these properties is the introduction of a suitable operational description. Except in 
the limit of infinite on-site Coulombian repulsion U 00 16, 1^ Il8l Il9l|. the introduction of such a description 
is an open problem of great physical interest. For low values of energy useful information about the effects of the 
Q non-perturbative electronic correlations is provided by two-component conformal- field theory [20I I21I |2^ I23L |2^ . 
Q Unfortunately, that method does not apply for finite values of energy. 

ILjl In view of the above-mentioned unusual finite-energy/frequency spectral properties observed in real experiments, 

which are far from being well understood, efforts towards the introduction of a suitable operational description to 
deal with the finite-ener gy pro blem are welcome. In this paper we introduce an operational representation for the ID 
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5h ' Hubbard model 25l I26ll27ll28| in terms of non-interacting pseudofermions. We find that in the thermodynamic limit 
. . 1 the wave function of excited states contained in the few-electron excitations factorizes into separated contributions 
corresponding to different pseudofermion branches. (A few-electron excitation is generated by application onto the 
ground state of operators whose expression involves the product of a few electronic creation and/or annihilation 
operators.) Such factorization occurs for all values of the on-site Coulombian repulsion U of the ground-state normal- 
ordered ID Hubbard model. The pseudofermion o per ational description is closely related to the pseudoparticle 
representation previously considered in the literature HE HE], and is the natural starting point for studies of the finite- 
energy /frequency few-electron spectral properties. As a result of the wave-function factorization, the pseudofermion 
description is more suitable for the study of the overlap between few-electron excitations and the energy eigenstates 
than the pseudoparticle representation. 

Our star ting point is a holon, spinon, and cO pseudoparticle representation, which refers to the whole Hilbert space of 
the model 31]. The relation between the original electrons and these elementary quantum objects involves the concept 
of rotated electron. The rotated electrons are related to the electrons by a unitary transformation first introduced in 
Refs. [sEll^- For such rotated electrons double occupation is a good quantum number for all values of U. Except 
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for the cO pseudoparticles, all pseudoparticle branches introduced in Ref. [29| have a composite character in terms 
of holons or spinons |3l| . The pseudofcrmions are related to the pseudoparticles by a second unitary transformation. 
The charge or spin carried by the pseudoparticle and its holon or spinon contents remains invariant under such a 
transformation. It follows that except for the cO pseudofermions, all remaining pseudofermion branches are composite 
objects of holons or spinons. The concepts of local pseudoparticle and effective lattice widely used in this paper are 
introduced in Ref. |35l |. 

The pseudofermion operational description introduced in this paper and the associated factorization of the wave 
function of energy eigenstates contained in few-electron excitations are used elsewhere in the study of the finite-energy 
spectral properties [iL ITol Is^. The theoretical predictions of Refs. P, ^3 seem to describe both qualitatively and 
quantitatively the one-electron removal spectral lines observed by photoemission experiments for finite values of the 
energy in real quasi-lD materials. We thus expect that our pseudofermion operational description is useful for the 
further understanding of the exotic properties displayed by low-dimensional materials. 

The paper is organized as follows: In Sec. II we introduce the ID Hubbard model and the rotated electrons. 
The elementary holon, spinon, and cO pseudoparticle description and associated au pseudoparticle representation are 
summarized in Sec. III. Moreover, in that section we introduce the operator algebra for the pseudoparticles and 
the statistics of these quantum objects according to the generalized Pauli principle In Sec. IV we consider the 
ground-state normal-ordered pseudoparticle operator description and introduce useful ground-state quantities. The 
pseudofermion description and the relationship between pseudoparticle and pseudofermion operators are introduced 
and discussed in Sec. V. This includes the introduction of the pseudofermion anticommutator algebra. In Sec. VI 
we study the pseudofermion energy and momentum spectra and introduce and discuss the factorization of the few- 
electron Hilbert subspace of the ground-state normal-ordered ID Hubbard model. The investigation of the laws under 
the pseudoparticle - pseudofermion transformation of several quantum objects and quantities is the subject of Sec. 
VII. Finally, in Sec. VIII we present the discussion and the concluding remarks. 



II. THE ID HUBBARD MODEL AND ROTATED ELECTRONS 



In a chemical potential fi and magnetic field H the ID Hubbard Hamiltonian can be written as, 



H = Hsoi4) + , (1) 

a—c, s 

where the Hamiltonian 



Hsoii) ^Hh- (C//2) N + (C//4) iV, ; Hh^T + UD, (2) 

has SO ( A) symmetry. Here Hh is the "simple" Hubbard model, 



T = -t 



(3) 



is the kinetic- energy operator, and 



Na Na 

^^Yl 4, T "^J' T 4- 1 ^i- i = Y T ri], i , (4) 

is the electron double-occupation operator. On the right-hand side of Eq. we have that /ic = fj., fis = fJ'oH, /io is 
the Bohr magneton, and the number operators, 



S!^--[Na-N]; 5f = --[iVT-7Vj, (5) 

are the diagonal generators of the ?7-spin and spin SU{2) algebras [s^ l3^ l39l|. respectively. We consider that the 
number of lattice sites Na is large and even and that Na/2 is odd. The electronic number operators on the right-hand 
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side of Eq. Q read N — J2cr=i i and Na- = J2j=i ^j,^^ where the operator Nj,a- = c] o-C^-, ^ counts the number 
of spin CT electrons at real-space lattice site j. On the right-hand side of Eqs. Q-© the operator ^ (and Cj^„) 
creates (and annihilates) a spin a electron at lattice site j — 1,2,..., Na- We denote the lattice constant by a and the 
lattice length hy L = Na a. 

The momentum operator is given by. 



^= E E^-(^)^=^ E r'^dkNa{k)k. (6) 

Here the spin a momentum distribution operator reads N^ik) = c\ ^Cfe^o-, where the operator ^ (and Ck,^) creates 
(and annihilates) a spin a electron at momentum k. The operators c\ ^ and Cfc_ ^ are related to the above operators 
^ and Cj^cr by the following Fourier transforms. 



Cfe.cr = 



(7) 



The Hamiltonian Hsn(A\ given in Eq. jSJ commutes with the six generators of the yy-spin and spin SU{2) algebras 
and has 5*0(4) symmetry {371 ISa . l39l|. While the expressions of the two corresponding diagonal generators are given 
in Eq. jSjl, the off-diagonal generators of these two SU{2) algebras read 



N, 



Sc 



E 



(8) 



and 



^l = E4 



Na 



(9) 



respectively. 

Throughout this paper we use units of Planck constant one and denote the electronic charge by — e. The Bethe- 
ansatz solvability of the ID Hubbard model © is restricted to the Hilbert subspace spanned by the lowest-weight 
states (LWSs) |23,|2^ or highest- weight states (HWSs) |23| of the 77-spin and spin algebras, that is by the states whose 
Sa and numbers are such that Sa = —S^ or 5'^ = 5^, respectively, where a = c for charge and a = s for spin. In 
this paper we choose the jy-spin and spin LWSs description of the Bethe-ansatz solution. In this case, that solution 
describes energy eigenstates with electronic densities n = N/L and spin densities m = [N^ — iVj/L in the domains 
< n < 1/a and < m < n, respectively. Some of our results correspond to the ranges < n < 1/a and < m < n. 
The description of the states corresponding to the extended domains 0<n<l/a ;l/a<n<2/a and —n < m < n 
; —{2/a — n) < m < (2/a — n), respectivel y, is achieved by application onto the latter states of off-diagonal generators 
of the 77-spin and spin SU{2) algebras |3lLl33|. 

Each lattice site j — 1,2, ...jNa of the model can either be doubly occupied, empty, or singly occupied by a 
spin-down or spin- up electron. The maximum number of electrons is 2Na and corresponds to density n — 2/a. Besides 
the N electrons, it is useful to consider [2Na — N] electronic holes. (Here we use the designation electronic hole instead 
of hole, in order to distinguish this type of hole from the pseudoparticle hole and pseudofermion hole.) Our definition 
of electronic hole is such that when a lattice site is empty, we say that it is occupied by two electronic holes. If a 
lattice site is singly occupied, we say that it is occupied by an electron and an electronic hole. If a lattice site is 
doubly occupied, it is unoccupied by electronic holes. The same definition holds for the rotated-electronic holes. We 
note that the lattice occupied by rotated electrons is identical to the original electronic lattice. 

The electron - rotated-electron unitary transformation maps the electrons onto rotated electrons such that rotated- 
electron double occupation, no occupation, and spin-up and spin-down single occupation are good quantum numbers 
for all values of U/t. We call cj ^ the electrons that occur in the ID Hubbard model © and while the operator c| ^ 
such that = V^U/t) ^ V{U/t) represents the rotated electrons, where the electron - rotated-electron unitary 
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operator V{U/t) is defined below. Similarly, ^ — ^i^/^)^] o-^^(t^/*)- Note that ^ and ^ are only identical 
in the U/t — > cxd limit where electron double occupation becomes a good quantum number. 

The operators V''{U/t) and V{U/t) associated with the electron - rotated-electron unitary transformation can be 
written as, 

V\U/t)^e-^; V{U/t)^e^. (fO) 
The operator S of Eq. pOfl is uniquely defined by the following two equations. 



Hh = V\U/t)HHV{U/t) = Hh + [Hh, S] + - [[Hh, 5], S] + ... , (ff) 

and 



[Hh, V^iU/t) D V{U/t)] - [Hh, D] ^ , (12) 
where the Hamiltonian Hh is given in Eq. ^ and the rotated-electron double occupation operator D reads. 



D^V^ {U/t) D V{U/t) = 4. T ^^-^ T 4, i c.-. i ■ (13) 

j 

Here D is the electron double occupation operator given in Eq. Q. The operator (|13|l commutes with the ID Hubbard 
model. We denote the rotated-electron double occupation by Dr- It is a good quantum number for all values of U/t. 

The transformation associated with the electron - rotated-electron unitary operator V{U/t) was introduced in Ref. 
[sH. The studies of that reference referred to large values of U/t and did not clarify for arbitrary values of U/t 
the relation of rotated-electron double occupation to the quantum numbers provided by the Bethe-ansatz solution. 
However, this transformation is uniquely defined for all values of U/t by Eqs. H1Q(I - (|12(I . Equations HlH and (|12|l can 
be used to derive an expression for the unitary operator order by order in t/U. The authors of Ref. (33| carried out 
this expansion up to eighth order (see foot note [12]). 



III. THE PSEUDOPARTICLE OPERATORS AND THE HOLON, SPINON, AND cO 
PSEUDOPARTICLE BASIC DESCRIPTION 



According to the studies of Ref. |23, there is an infinite number of pseudoparticle branches: the cO pseudoparticles 
and the av pseudoparticles such that a = c, s and ly = 1,2, .... The aiy pseudoparticle notation considered in this 
paper is related to the c pseudoparticle and 01,7 pseudoparticle notation of Ref. as follows: cO = c for = 0, 
J/ = 7 and cv = c,"f for v = 1,2, and = 7 -f 1 and sv = s,"/ + 1 for v = 1,2, .... Moreover, we denote by bare 
momentum q the momentum carried by the pseudoparticles, while in Ref. p9| it was called band momentum. Our 
designation is justified by the form of the pseudofermion momentum, as further discussed below. Note that within 
our notation, the general designation of au pseudoparticle refers to the a — c branches such that = 0, 1, 2, ... and 
a — s branches such that v = 1, 2, .... Elsewhere it is shown that the cv pseudoparticles and pseudoparticles are 
2j/-holon and 2i^-spinon composite objects, respectively, where = 1, 2, ... [sj- 

The introduction of the pseudofermion operational description studied in this paper requires the use of an operator 
representation for the pseudoparticles. Thus, in this section we introduce an operational description for the av pseu- 
doparticles and the holons and spinous which are not part of composite pseudoparticles. An operational description 
for the pseudoparticles of bare-momentum q was introduced in Ref. |2^. However, such a description did not take into 
account the holon (and spinon) composite character of the cv pseudoparticles (and sv pseudoparticles). Moreover, the 
spinon description used in the studies of that reference is not valid for the whole Hilbert space. This affected the values 
of the entries of the pseudoparticle statistical-interaction matrix j36l | . In this section we provide the correct values for 
the entries of such a matrix. Another limitation was the lack of a representation for the pseudoparticle operators in 
terms of spatial coordinates. The concepts of local av pseudoparticle and effective av lattice are introduced in Ref. 
[35I I . The studies of that reference include a description of the local pseudoparticles in terms of the rotated-electron 
site distribution configurations. The absence of such a local operational description is one of the reasons why the 
studies of Ref. |2^ did not provide useful information about the main issue of the relation between the pseudoparticle 
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and electronic operators. Also the concept of rotated electron is important for the study of that issue and was not 
considered in the studies of such a reference. All these concepts are valuable and necessary for the application of the 
operational pseudofcrmion representation introduced in this paper to the evaluation of few-electron spectral functions 

In Appendix A we summarize the basic properties of the av pseudoparticles which are needed for our studies. 
This includes introduction to the av pseudoparticle bare momentum, effective lattice, and useful ground-state 
quantities. 



A. THE PSEUDOPARTICLE OPERATORS AND PSEUDOPARTICLE STATISTICS ACCORDING TO 

THE GENERALIZED PAULI PRINCIPLE 

Generation and removal of pseudoparticles is in general associated with creation and/or annihilation of electrons. Yet 
there are also transitions which change the numbers of these quantum objects at constant spin a electron numbers. One 
can introduce elementary operators for creation or annihilation of av pseudoparticles. In this subsection, we introduce 
two alternative representations corresponding to pseudoparticle operators, both in terms of the bare-momentum q 
and spatial coordinates. These representations refer to the bare- momentum pseudoparticles and local pseudoparticles 
[ssf . respectively. 

Let us introduce the bare- momentum av pseudoparticle creation (and annihilation) operator 6 J ^.^ (and hq^av) 
which creates (and annihilates) a av pseudoparticle of bare momentum q. In addition, we introduce the local av 
pseudoparticle creation operator b\.. and annihilation operator bx^ , au ■ These bare-momentum and local pseudopar- 
ticle operators are related as follows. 



j = l j = l 

The local av pseudoparticle creation (and annihilation) operator h\.. (and bxj,av) creates (and annihilates) a local 
av pseudoparticle at the effective av lattice site of spatial coordinate Xj = aau j, where j — 1,2,..., N*^^. The effective 
av lattice constant is defined in Eq. (|A7I) of Appendix A, where the concept of an effective av lattice, introduced 
in Ref. |35j| is described. The conjugate variable of the bare-momentum qj of the av pseudoparticle branch is the 
space coordinate Xj of the corresponding effective av lattice. This is different to the electronic operators of Eq. |(7J), 
where the conjugate variable of the momentum kj is the space variable of the original electronic lattice. In reference 
|35j| . the pseudoparticle site distribution configurations in the effective av lattices are related to the corresponding 
rotated-electron site distribution configurations. 

The pseudoparticles obey a Pauli principle relative to the bare- momentum occupancy configurations, i.e. a discrete 
bare-momentum value qj can either be unoccupied or singly occupied by a pseudoparticle. For v > composite 
av pseudoparticles, the number of discrete momentum values N*^ is not the same for all energy eigenstates. Thus, 
these objects cannot be classified as fermions or bosons. In order to classify the av pseudoparticle according to the 
generalized Pauli principle introduced in Ref. |36j . we consider the av dimensions, 



= 1 + N*, - 7V„. = 1 + N^, , (15) 

where according to Eqs. (|A3p - (|A6|) of Appendix A, Nai, and A^^^^ are the number of av pseudoparticles and av 
pseudoparticle holes, respectively, and N*^^ = Nai, + N^^^. A transition to an excited energy eigenstate produc- 
ing deviations l^Na^ in the av pseudoparticle numbers, leads to the following deviations in the corresponding av 
dimension. 



Adajy = - ^ ^ gav,a'v' ^Na'u' ■ (16) 
q'— c, s u'—6^i ^ 

According to the generalized Pauli principle 36], the parameters gav,a'y' on the right-hand side of this equation are 
the entries of the statistical-interaction matrix. From the use of Eqs. ljA3|) - (IA6|) of Appendix A, we find that for the 
av pseudoparticles such a statistical-interaction matrix has infinite dimension and its entries are given by, 

gcO,av = 5a,c5v,0\ Qav, cO = ^a, c - 5a, s Sai/, a V = (5^, ( + I^' - 1 1/ - J/' | ) ; V,v'>Q. (17) 
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This fully defines the statistics of the au pseudoparticles. We emphasize that the cO pseudoparticle entry Qco.av given 
in Eq. IjlTI) has Fermionic character. This is related to the fact that the number of sites of the effective cO lattice 
is constant and given by Na, and that the corresponding effective lattice constant Oco equals the electronic lattice 
constant a. Thus, the effective cO lattice and electronic lattice are identical. Furthermore, it is found in Ref. ^| 
that the sites occupied by cO pseudoparticles (and cO pseudoparticle holes) are the same as the sites singly occupied 
by rotated electrons (and doubly occupied and unoccupied by rotated electrons). While the cO pseudoparticles do 
not feel the statistical-interactions of the remaining pseudoparticles, the v > Q composite av pseudoparticles feel the 
statistical interactions of the cO pseudoparticles, as confirmed by Eq. ((T7|) . The form of the entries given in Eq. l(T7|) 
confirms that the latter composite pseudoparticles are neither fermions nor bosons "s^. The statistical interactions 
of the composite pseudoparticles result in part from the property that the width of the bare-momentum domain 
SfZaiy — 27r[l/aQ,y — 1/ L] is different for different values of the pseudoparticle numbers. 

The au pseudoparticle bare- momentum distribution functions Nau{q) play an important role in the pseudoparticle 
description (soi, Isij . These functions are for all energy eigenstates the eigenvalues of the following pseudoparticle 
bare-momentum distribution operators. 



NcM^blc^uKo..- (18) 
The bare-momentum distribution functions Na_v{q) read Nau{<lj) ~ 1 for occupied discrete bare- momentum values qj 
and Na,y{qj) ~ for unoccupied discrete bare- momentum values qj. Each LWS is uniquely specified by the values 
of the set of distribution functions {Naiy{q)} such that v — 0, 1,2, ... for a — c and ly — 1,2, ... for a = s. Physical 
quantities such as the energy, depend on the values of these distribution functions and numbers through the rapidity 
momentum functional k{q) and rapidity functionals Acy(g') and Ac^oV The value of these functionals is uniquely 
provided by solution of the following functional integral equations |3Ct l3l| , 



kcu{q) = q-\ — / dq A'co(9 ) arctan' 



vU/At 

iY.r ■'^'"-(■'^^..(i^iiiy^); .>0, (20) 



and 



n 1 ^ 'AT f , ^A..(g)-sinfc(g') 

= q / dq iVco(g ) arctan' 



iyU/4:t 



^E£>^-^.^^tf)e..( ^-'%:-'^'' ). (21, 



Here 



kc^iq) = 2Re{arcsin(Ac^((7) + ivU/U)} ; > , (22) 

is the cv rapidity-momentum functional and the limiting bare-momentum values 5^ and qau where a = c, s and 
1/ = 1,2, ... are given in Eqs. (|A9|I and HA10|I and in Eq. (|A11() . respectively, of Appendix A. The function Q^^^i{x) 
appearing in Eqs. H2U|I and (|21|l is given in Eq. (|B17|) of Appendix B. The equations (|19|I - H21II correspond to a 
functional representation of the thermodynamic Bethe-ansatz equations introduced by Takahashi [26j . The rapidity- 
momentum functional is real and the rapidity functionals are the real part of Takahashi's ideal strings |2a. IsH . It is 
useful to introduce the following cO rapidity functional. 



Aco(q) = sinfc(g) , 

where k{q) is the rapidity-momentum functional. 



(23) 
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B. HOLONS, SPINONS, cO PSEUDOPARTICLES AND THE YANG HOLON AND HL SPINON 

ELEMENTARY OPERATORS 

In this subsection we consider a holon, spinon, and cO pseudoparticle description for the whole Hilbert space of the 
ID Hubbard model whose validity is shown elsewhere _31J. Such elementary quantum objects correspond to specific 
rotated-electron site occupations. This description of all energy eigenstates in terms of occupancy configurations of 
three elementary quantum objects only, is useful for the studies of this paper. 

Let us distinguish the total rj spin (and spin) value, which we denote by Sc (and Sg) and the corresponding 77-spin 
(and spin) projection, which we denote by (and S^), from the rj spin (and spin) carried by the elementary quantum 
objects. We call Sc (and Ss) the rj spin (and spin) carried by the holons, spinous, and other elementary objects and dc 
(and as) their rj-spm (and spin) projection. The operators Mc,a^ and Ms, o-, which count the number of the CTc = ±1/2 
holons and as = ±1/2 spinons have the following form [sij . 

Mc,-i/2 = VHU/t)Y,<^hc,Tc]^Cj^ViU/t); 
j 

Me, +1/2 = V\U/t) J2 c.T 4t V{U/t) , (24) 

j 

and 



Ms,-i/2 = VHU/t)J2clc,^c]^c,^V{U/t); 
j 

Ms, +1/2 = VHU/t) J2 4t ^J-T ^(U/t) , (25) 

j 

respectively. Here, the operator c^-|. Cj| cj^ Cji counts the number of electron doubly-occupied sites, 

''■'T "^]t '^^^ counts the number of electron empty sites, cj^ Cj| cj| Cji counts the number of spin-down electron 

singly-occupied sites, and c'j^ cji c^j^ Cj| counts the number of spin-up electron singly-occupied sites. The operator 

V{U/t) on the right-hand side of Eqs. and ((53 is uniquely defined for all values of U/t by Eqs. (^-(^3). The 
new physics brought about by the relations of Eqs. H24|l and (|25|l is that the electron - rotated-electron unitary trans- 
formation generates the holons and spinons whose occupancy configurations describe the exact energy eigenstates. 
Throughout this paper we denote such holons and spinons according to their value of ac — ±1/2 and as = ±1/2, 
respectively. The unitary rotation is such that for all values of U /t the number of emerging —1/2 holons, +1/2 holons, 
— 1/2 spinons, and +1/2 spinons equals precisely the number of rotated-electron doubly occupied sites, empty sites, 
spin-down singly occupied sites, and spin-up singly occupied sites, respectively. 

The holons have Sc = 1/2, Ss = 0, and as = 0. Thus, there are two types of holons, which have CTc = ^1/2 and 
CTc = +1/2 and carry charge — 2e and +2e, respectively. These objects are on-site spin-singlet rotated-electron pairs 
and rotated-electron unoccupied sites, respectively. An important property is that rotated-electron double occupation 
Dr equals the value M^^_i/2 for the number of —1/2 holons. The rotated-electron double occupation Dr = M^,, _i/2 
plays an important role in the description of the few-electron spectral properties of the quantum problem. For few- 
electron spectral functions the weight distribution resulting from transitions to the Hilbert subspace spanned by excited 
states of rotated-electron double occupation Dr corresponds to the Dl!^ upper Hubbard band, where Dr = 1,2,... 
|34j . Moreover, the spin-up and spin-down rotated-electron singly occupied sites correspond to the as = +1/2 and 
as = —1/2 spinons, respectively. The Sg — 1/2 spinons have no charge degrees of freedom and thus describe the 
spin degrees of freedom of rotated-electron singly occupied sites only. The charge degrees of freedom of the rotated 
electrons (and rotated-electronic holes) of these singly occupied sites are described by the chargeons of charge — e 
(and antichargeons of charge +e) [sj- These quantum objects are part of the charge degrees of freedom, but do 
not contribute to the 77-spin and spin SU{2) algebras. In the case of the description of the transport of charge in 
terms of electrons (and electronic holes) the elementary carriers of charge are the chargeons of charge — e and —1/2 
holons of charge — 2e (and the antichargeons of charge +e and +1/2 holons of charge +2e). The cO pseudoparticle 
has no ry-spin and no spin degrees of freedom and is a composite quantum object which contains a chargeon and 
an antichargeon. However, in how transport of charge is concerned, the chargeon and antichargeon correspond to 
alternative descriptions of the cO pseudoparticle. When the transport of charge is described in terms of electrons (and 
electronic holes) the cO pseudoparticle couples to charge probes through the chargeon (and antichargeon) and carries 
elementary charge — e (and +e). 
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All energy eigenstates of the ID Hubbard model can be described in terms of occupancy configurations of holons, 
spinous, and cO pseudoparticles. On the other hand, in Ref. it was found that all energy eigenstates associated 
with the ID Hubbard model Bethe-ansatz solution |2Ml2a|' be described in terms of occupancy configurations ofav 
pseudoparticles. By merging these two representations, one finds that the cv pseudoparticles such that i/ ^ and the 
sv pseudoparticles are 2v- holon and 2.-spinon composite objects, respectively M- The composite c. pseudoparticles 
(and sv pseudoparticles) are Sc ^ Uc — (and Sg = Cs = 0) objects without spin (and charge) degrees of freedom. 
These composite quantum objects contain an equal number v of —1/2 holons and +1/2 holons (and —1/2 spinous and 
+1/2 spinous). When the transport of charge is described in terms of electrons (and electronic holes), the 2i/- holon 
composite pseudoparticles couple to charge probes through their ly —1/2 holons of charge — 2e (and their +1/2 
holons of charge +2e). Thus these composite quantum objects carry elementary charge — 2i/e (and +2iye). 

The ±1/2 holons (and ±1/2 spinous) which are not part of 2i^-holon composite pseudoparticles (and 2i/-spinon 
composite si/ pseudoparticles) are called ±1/2 Yang holons (and ±1/2 HL spinous). In the designations HL spinon 
and Yang holon, HL stands for Heilmann and Lieb and Yang refers to C. N. Yang, respectively, who are the authors of 
Refs. 01 Ell- Note that the holons (and spinous) which are part of Sc = (and Ss = 0) composite a/ pseudoparticles 
(and siy pseudoparticles) remain invariant under application of the off-diagonal generators of the ry-spin (and spin) 
SU{2) algebras given in Eq. (and l|njl). Application of these generators produces 77-spin (and spin) flips in ±1/2 
Yang holons (and ±1/2 HL spinous). Thus, the —1/2 Yang holons and —1/2 HL spinous are created by application 
onto the LWSs of suitable off-diagonal generators of the ry-spin and spin SU{2) algebras, respectively, given in Eqs. 
© and The generators © and can be written in terms of rotated-electron operators as. 



-1/2 



iV, 

E 



-1/2 



(26) 



and 



Na Na 

respectively. In these equations qc ~ n/a and = is the momentum. The operators (|26f) and (|27|l are invariant 
under the electron - rotated-electron unitary transformation j3l| |. Therefore, they have precisely the same expression 
in terms of electron and rotated-electron creation and annihilation operators. 

Note that within the operational representation of Eqs. H26II and H27() the Yang +1/2 holons and HL +1/2 spinons 
are not explicitly considered. Indeed, the numbers ic, +1/2 of +1/2 Yang holons and is.+i/2 of +1/2 HL spinons 
are fully determined by the numbers ^c, -1/2 of Yang —1/2 holons, is,-i/2 of HL —1/2 spinons, and the set {Nav} 
of the different av pseudofermion branches. This justifies why here we consider the operators of Eqs. (|26|) and H27|) 
as creation and annihilation operators for —1/2 Yang holons and —1/2 HL spinons, respectively. When applied onto 
LWSs, the operators -^^2 produce energy eigenstates with finite values for the numbers ^c, -1/2 and/or is,-i/2 
of the following form. 



|{i,,_i/2,i.,-i/2})- n '^'--"tS' ^'^ LWS). (28) 

a — c. s V 

Here \LWS) is the LWS that corresponds to the state |{-Z^c,-i/27 ^s, -1/2}) and La ~ L^, +1/2 + La, -1/2 — "^Sa where 
Sa is the state T^-spin value (a = c) and spin value (a = s). The energy eigenstates (|28ll are not described by the 
Bethe-ansatz solution [s^ • 

The pseudoparticle bare momentum q is such that q G (—qav, +qaiy) where the limiting bare-momentum values iqav 
associated with the limits of the ai/ pseudoparticle Brillouin zone are given in Eq. IjAlip of Appendix A. According 
to the results of Ref. j3fl], the ±1/2 Yang holons, ±1/2 HL spinons, and v > Q composite av pseudoparticles of 
bare-momentum values q = ±qai> are non-interacting and localized quantum objects. Such behavior results from the 
invariance of these quantum objects under the electron - rotated-electron unitary transformation. This means that the 
±1/2 Yang holons, ±1/2 HL spinons, and v > composite av pseudoparticles of bare- momentum values q — iqav 
are the same quantum objects as the corresponding rotated objects. Thus, these objects are localized and do not 
contribute to the transport of charge or spin. However, in general a v > composite av pseudoparticle is different 
from a, V > rotated composite av pseudoparticle. The only exception is precisely for bare momentum values q 
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such that q -i^qav As the bare momentum approaches its limiting values, q iqai^, the u > composite ap 
pseudoparticle and the v > rotated composite ai' pseudoparticle become the same quantum object. Importantly, 
the electrons and rotated electrons involved in processes associated with creation of ±1/2 Yang holons, ±1/2 HL 
spinous, and > composite av pseudoparticles of bare-momentum values q = ±qaD are also the same quantum 
object, i.e. remain also invariant under the electron - rotated-electron unitary transformation. This always refers to 
localized electrons. It follows that the transport of charge (and spin) is associated with the cO pseudoparticle and 
q ±qcv composite pseudoparticle quantum charge fluids (and q ^ ±qsi> composite sv pseudoparticle quantum 
spin fluids) where v — 1,2, .... 

Within the pseudoparticle, Yang holon, and HL spinon operational description, the ±1/2 holon (a = c) and ±1/2 
spinon (a = s) number operators M^^±i^2 given in Eqs. (|24|l and H25|l are written in terms of pseudoparticle operators 
and ±1/2 Yang holon (a — c) or ±1/2 HL spinon (a — s) number operators ^^,±1/2 as follows, 



00 +qa,u 

M„,±i/2 =L„,±i/2+^ i^Nc^.iq), (29) 

where the pseudoparticle bare-momentum distribution operators Nau{q) are provided in Eq. (|18|) . The operator 
La, ±1/2 can be written as, 



La, ±1/2 = V ^c..L + 1/4 - 1/2 T SI ; a = c,s. (30) 

Here is the diagonal generator of the 77-spin (a — c) and spin (a — s) algebras whose expression in terms of 
electronic operators is provided in Eq. ((S)). The spin a electron number operator commutes with the electron - 
rotated-electron unitary operator. Thus, the operator has the same expression in terms of electron and rotated- 
electron creation and annihilation operators. The same occurs for the ?7-spin {a — c) and spin (a = s) operator 

Sa.Sa. This operator can be expressed in terms of rotated-electron creation and annihilation operators by use of Eqs. 
(|26|l and (|27l) . Therefore, Eq. (|3U|I provides the expression of the operator La,±i/2 both in terms of electron and 
rotated-electron creation and annihilation operators. 

An electronic ensemble space is spanned by all energy eigenstates with the same values for the numbers Nf and 
Ni . An important concept is that of CPHS ensemble space 30] . This is a subspace spanned by all energy eigenstates 
with the same values for the numbers of ±1/2 holons and ±1/2 spinous {Ma j-1/2}, where a = c,s. In general, an 
electronic ensemble space contains several CPHS ensemble spaces. Moreover, usually a CPHS ensemble space includes 
different CPHS ensemble subspaces. A CPHS ensemble subspace is spanned by all energy eigenstates with the same 
values for the sets of numbers Ncq, {Naiy}, and {^0,-1/2} such that a — c,s and i' — 1,2,.... (According to the 
notation of Ref. [sOl, CPHS stands for c pseudoparticle, holon, and spinon.) 



C. THE GENERAL ENERGY SPECTRUM AND THE MOMENTUM OPERATOR 



The generators of the ry-spin and spin SU{2) algebras given in Eqs. (|26|l and H27|l commute with the cO pseudoparticle 
and composite av pseudoparticle creation and annihilation operators of Eq. (I14|l . We emphasize that this property 
has the following important effect: All 2Sa energy eigenstates obtained from a given regular energy eigenstate have the 
same pseudoparticle momentum distribution functions Nca{q) and {Nauiq)} for all branches a = c, s and = 1, 2, .... 
Thus these states are described by similar pseudoparticle occupancy configurations and only differ in the relative 
numbers of +1/2 Yang holons and —1/2 Yang holons {a = c), or/and in the relative numbers of +1/2 HL spinous 
and —1/2 HL spinous (a = s). This reveals that the coupled functional equations (|19|l - (|21|) . which involve the 
pseudoparticle momentum distribution functions and do not depend on the La, numbers, describe both LWSs and 
non-LWSs. The eigenvalues E of the energy eigenstates of the Hamiltonian |^ can be written in the following form. 



a—c, s 

where 



E 



SO(4) 



Mc-2M,,_i/2-^ 



(32) 



10 



and 



Eh = -2t— / dqNcoiq) cos k{q) 



OO i-+q, 



At—Y 



dq N,,{q) Re { - (Ac.(g) + U/Atf'^ + U L,, _i/2 . (33) 



On the right-hand side of Eq. the numbers = —\[Na — N] and = — ^[^t — N^] are the eigenvalues of the 
diagonal generators given in Eq. 101 and /ic = 2/i and fis — 2/io H are the same quantities as on the right-hand side 
ofEq. CEJ. 

The values of the rapidity- momentum functional k{q) and rapidity functionals Aai^iq) are the same for all 25*^ -I- 1 
states in the same tower. Such functionals are eigenvalues of operators which commute with the off-diagonal generators 
of 77-spin and spin algebras. This is consistent with the ry-spin and spin SU{2) symmetries, which imply that the 
Hamiltonian ^ commutes with these generators and thus the energy (|32|l is the same for the set of 2Sa + 1 states 
belonging to the same 77-spin {a = c) or spin {a ~ s) tower. The above operators obey equations similar to Eqs. H19() - 
()21|l . with the av bare-momentum distribution functions Nav{q) replaced by the corresponding operators Nav{q) given 
in Eq. H18(l . The pseudoparticle Hamiltonian expression is also obtained by replacing in the energy expressions H31(l - 
(I33|l the distribution functions Nav{<i} by the operators Nav{q)i and the rapidity- momentum and rapidity functionals 
by the corresponding operators. 

The momentum operator given in Eq. Q can be expressed in terms of pseudoparticle and —1/2 holon operators 
as follows ||31j . 



9^0 00 +qsv 00 +9c 



P= Y.N,oiq)q + Y^ J2 N,,{q)q + J2 E NUl) ~ l] + ^ -1/2 

q=q-^ y=lq=-qs„ l'=lq=-qa,„ 

qto 00 +q3v 00 u 

- E ^co(<z)g + E E ^-(9)9 + E E ^c.(<z) [(1 + ^^)^-91 + ^4-1/2- (34) 

q=q-^ u=lq=-qsv u=lq=-qc,u 

The momentum operator (|34|l commutes with the electron - rotated-electron unitary operator V{U /t) j31|. The mo- 
mentum eigenvalues can be straightforwardly written by replacing on the right-hand side of Eq. H34(l the pseudoparticle 
and holon number operators by the corresponding eigenvalues. 

Note that the number -1/2 and corresponding operator ic, -1/2 appearing on the right-hand side of Eqs. H33|l 
and H34|l . respectively, are beyond the Bethe-ansatz solution. Indeed, the states which span the Hilbert subspace 
associated with that solution have no —1/2 Yang holons and thus Lc, -1/2 = in Eq. (|33|l for these states. However, 
the energy spectrum defined by Eqs. (|31|l - l|33|) and the momentum-operator expression (|34l) refer to the whole Hilbert 
space of the ID Hubbard model. On the other hand, as a result of the SO{A) symmetry of the Hamiltonian the 
corresponding energy spectrum (|32|) docs not depend on the value of the number L^, -1/2- This property is confirmed 
by analysis of expressions and and by noting that Mc_„i/2 = ^c, -1/2 + Sl^i ^^cy 



IV. THE GROUND-STATE NORMAL-ORDERED PSEUDOPARTICLE OPERATIONAL 

DESCRIPTION 

In this section we consider the ground-state normal-ordered pseudoparticle operational description which is needed 
for the introduction of the pseudofermion description. We are mostly interested in excited states generated from the 
ground state by processes involving changes in the occupancy configurations of a finite number of av pseudoparti- 
cles, —1/2 Yang holons, and —1/2 HL spinous. The pseudofermion description introduced in the ensuing section 
corresponds to the Hilbert subspace spanned by such excited states. In the thermodynamic limit the few-electron 
excitations are contained in that subspace. 



A. NORMAL-ORDERED NUMBER OPERATORS AND USEFUL GROUND-STATE QUANTITIES 



Throughout this paper the symbol : O : refers to the ground-state normal-ordered expression of a general operator 
O. Such a ground-state normal-ordered expression is given by that operator minus its ground-state expectation value. 
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We introduce the ground-state normal-ordered av pseudoparticle bare-momentum distribution operator, 

■NaAq):^NaM-KAi). (35) 
and the —1/2 Yang holon (a = c) and —1/2 HL spinon (a = s) ground-state normal-ordered number operator, 

'■ La, -1/2 '■— La, -1/2 — L^a, -1/2 = La, -1/2 • (36) 

Here the operators Nau{<l) and La, -1/2 are defined in Eqs. (|18|1 and J^OJ), respectively, N^^{q) is the ground-state 
pseudoparticle distribution function whose expressions are given in Eqs. (jA12|) - (|A14|) of Appendix A, and L° — 
is the —1/2 Yang holon (a = c) and —1/2 HL spinon {a = s) ground-state number. It follows from the expressions 
given in Eq. (jA15(l of Appendix A for the ground-state values of the number iV*^, whose general expressions are 
given in Eqs. ljA3|l . (|A4|) . and (IA6p of the same Appendix, that in the case of the ground state the effective aiy lattice 
constants (|A7|) are given by ;35j. 



'^cO — ^ i '^Cf — sr i '^sl — J '^Sl/ — J (37) 

rif m 

where 6 — (1/a — n) is the doping concentration. We note that the meaning of the divergences in the value of the 
constants defined in Eq. l|57|l is that the corresponding effective av lattice has no sites, i.e. N^'* = and, 
therefore, does not exist for the ground state. This is the case of the effective cv lattices for half filling when v > 
and of the effective sv lattices for zero spin density when v > 1. It follows that such singularities just indicate the 
collapse of the corresponding effective av lattice. This is one of the reasons why some of our expressions refer to 
electronic densities < n < 1/a and spin densities < m < n, such that all ground-state effective av lattice constants 
(|37|l have finite values. 

From use of expressions (|S7|l one can write the ground-state number NaC given in Eq. HA15|) of Appendix A as, 



A^r-i-- (38) 

Except for terms of order 1/L, the limiting bare-momentum values defined in Eqs. ljA8p - ljAll|) of Appendix A simplify 
and are given by, 



(39) 



where the ground-state effective-lattice constants a^^ are these given in Eq. (|37|l . Comparison of Eqs. (|37|l and H39|l 
leads to the following ground-state expressions for the av pseudoparticle limiting bare-momentum values. 



= [7r/a-2kF], v>0; ql ^ [kp^ - kpi] , v>l. (40) 

In most situations, one can disregard the l/L corrections and use the bare-momentum limiting values given in Eq. 
(gOJ. 

For each specific energy eigenstate the rapidity- momentum functional k{q) and the rapidity functionals Aau{q) 
become mere functions of q and are called rapidity-momentum and rapidity functions, respectively. We denote the 
ground-state rapidity- momentum function by k°{q) and we call A'^i,{q) and A°j,(g) the ground-state rapidity functions. 
These ground-state functions are computed by solution of the integral equations obtained by introducing in Eqs. H19|) - 
(|21ll the ground-state distribution functions (|A12|) - (jA14p given in Appendix A. The solution of these equations can be 
written in terms of the inverse functions of k'^{q) and A^j^(q) which we call qc{k) and qau{A.), respectively. Here k and 
A are the rapidity-momentum coordinate and rapidity coordinate, respectively. For the pseudoparticle description 
these functions refer to the ground state only. However, in the case of the pseudofermion description introduced in 
the ensuing section these functions apply both to the ground state and excited states. It follows that these functions 
play an important role in the latter description. These are odd functions such that, 



qc{k) = -qd-k); qau{A) ^ -q^^{~A) ; a ^ c, s ; v ^ 1,2, 



(41) 
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The domains of the rapidity- momentum coordinate k and rapidity coordinate A are such that —it /a < k < +7r/a and 
—00 < A < oo, respectively. If follows that, 

fco(±^)=±^; A«o(±^)=0; AL(±C)=±oo; a = c, 5 ; i. = l,2,..., (42) 
where we introduced the ground-state cO rapidity function, 

A^o(<z) = sinfc"(q). (43) 

The cO rapidity function H43() is the ground-state value of the rapidity functional defined in Eq. (|23|l . The relations 
(|^ are equivalent to, 

gc(±^)=±^; gc.(±oo) =±[^-2fcF]; gs.(±(») =±[fcfT-fcFi]; j^ = 1,2,.... (44) 

The Fermi points ±2fci? and ±kpi of the cO pseudoparticles and si pseudoparticles, respectively, correspond to the 
following rapidity Fermi points, 

Q = k\2kF), B = AUkFi). (45) 

Then, 

qci±2kF) ^ ±Q ; qsi{±kFi) = ±B. (46) 

Interestingly, one can write the above functions qdk), qci/(A), and gs,y(A) in terms of the two-pseudofcrmion phase 
shifts expressed in terms of the rapidity-momentum coordinate k and rapidity coordinate A, whose physical 

meaning is discussed in the ensuing section. These functions read, 

Uk)^k + / &'$,o,co(fc',A:) , (47) 
J-Q 



qcu{A) = 2Re{eiicsm(A + ivU/4t^} - J dk' ^c.o,cuik' ,K) ] i/ = l,2,..., (48) 



f+Q 

qsAA)= dfc'$,o,s.(fc',A) ; i/ = l,2,..., (49) 
■J-Q 

respectively. The two-pseudofermion phase shifts play an important role in the few-electron spectral properties and 
appear in the expression of the momentum carried by the pseudofermions. The two-pseudofermion phase shifts 
^af, a'v' expressed in terms of the rapidity-momentum coordinate k and rapidity coordinate A are mathematically 
defined by the following equations, 

- /4isinfc 4tsinfc'\ ~ , - / U sink 4tA'\ , , 

$cO,co(fc,fc ) = *cO,cO ( ^ , ^ 1; «'cO,a;.(fc,A') = $c0.a^. ( ^ ' "[T j ' ^^^^ 

/X - /4iA 4tsinfc'\ ~ , , . - f^tA AtA'\ , , 

5>ai.,c0(A, fc') = $a<.,cO ( — , JJ j; ^ au, av' {A, A') = ^ av., av' \— , \ , (51) 

where the two-pseudofermion phase shifts ^av, a'v' {t, '"') are defined by the integral equations ljB2p - (|B15|) of Appendix 
B. These are expressed in terms of the variable r defined in the arguments of the functions on the right-hand side of 
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Eqs. I|5U|I - I|51|) . The values of the parameters Q and B introduced m Eq. I|45|) are controlled by the two-pseudofermion 
phase shifts and are computed by self-consistent solution of the following equations, 

2kF^Q+ dk<^>co,co{k,Q) ; kFi= dk <^>co, si {k, B) , (52) 

J-Q J-Q 

respectively. At zero spin density the parameter Q changes from Q ^ kp in the limit U/t Q to Q — 2kp as J7 — > oo. 
At half filling and zero spin density, the U/t dependence of the parameter Q is singular at U/t = 0. It reads Q — -k /2a 
at U/t — and is given by Q = "'/'J for all finite values of U/t. The U/t > value, Q = 2kp = n/a, is associated 
with a full cO pseudoparticle band both in bare-momentum q and in rapidity-momentum k spaces and thus implies 
insulator behavior. In contrast, the U/t = value Q = 7r/2a corresponds to a metallic band which is half filled in 
the rapidity-momentum space. Such a singular behavior is associated with the Mott-Hubbard transition 0,113. At 
zero spin density the parameter B is given by _B cx3 and vanishes in the limit of spin density to ^ n, as the fully 
polarized ferromagnetic state is approached. 

Unfortunately, it is in general difficult to obtain closed form expressions for the ground-state functions k'^{q) and 
^01/(9) by inverting the functions defined in Eqs. (|47|) - H49|I . This can be achieved for specific limits of the parameter 
space only. For instance, from the use of Eqs. (|47|) - (|49|l we can obtain for zero spin density to = 0, values of the 
electronic density < n < 1/a, and limiting on-site repulsion values U/t ^ and U/t >> 1 the following closed form 
expressions for the ground-state functions k^{q), A°g((7), A^^{q), and A^-^^{q), 



k'^iq) = ^; \q\<2kp, U/t-^O: 



= sgn(g) - kp] 
= sgn(g)7r/a; 



2kp < \q\ < Tr/a , 
\q\=7r/a, U/t^O; 



— ln(2) ^^^^^ 
U a 



U/t^O; 
q\ < Tr/a , U/t » 1 , 



(53) 



A2o(g) = sin(^); \q\ < 2kp , U/t^O; 

= sgniq) sm(^{\q\ ~ kp) 2kp < \q\ < n/a , U/t^O; 

= 0; \q\=7r/a, U/t^O; 

2tn 

= sin(ga) - — ln(2) sin(2ga) ; \q\ < n/a, U/t » 1 , (54) 



Aliq) ^ sgn(<z) sin(M_l_^) ; < |g| < (^/a - 2kp) , U/t^O 

= 0; g==0, U/t^O 

= ±00; q^±{Tr/a-2kp), U/t^O 

= ^tan(||); < \q\ < [n/a - 2kp) , U/t » 1 , (55) 

for h' > 0, and 



AOi(g) = sin(ga); \q\ < kp , U/t^O 
= ±00; q^±kp, U/t~^0 

= -^arcsinhftan( — )V \q\ < kp , U/t»l, (56) 
nU V 71 / 

respectively. We note that for zero-spin density the v > \ ground-state rapidity function A^^, (g) vanishes for all values 
of n and U/t. For to = 0, the bare-momentum domain width of the sv pseudoparticles belonging to branches such 
that V > \ vanishes and the corresponding effective sv lattice has no sites and thus collapses. This collapsing is 
associated with the divergence of the effective sv lattice constant given in Eq. (|37|l as to ^ when v > 1. The same 
occurs for the cv pseudoparticles belonging to branches such that v > when half filling is approached. 
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For any energy eigenstate {ip) the normal-ordered bare- momentum distribution and number operators defined in 
Eqs. l|S5|l and respectively, obey the following eigenvalue equations, 

: N^,{q) : lij) ^ AN^,{q)\,p) ; : L^, _,/2 : W = AL^, . (57) 

Here ANa^(q) is the pseudoparticle bare-momentum distribution function deviation and AL^ ~i/2 is the deviation 
in the number of —1/2 Yang holons (a = c) or of —1/2 HL spinous (a — s). These deviations are given by 

AiV„,(g)=iV„,(g)-iV°,(g); AL„, _i/2 = L., -1/2 - i° , -1/2 ■ (58) 

These values describe deviations of occupancy configurations of excited states relative to the ground-state occupancy 
configurations described by the bare- momentum distribution functions and numbers given in Eqs. (|A12p - l|A14|l of 
Appendix A. For these excited states, the at/ pseudoparticle bare-momentum distribution function and the —1/2 
Yang holon and —1/2 HL spinon numbers read, 

7Va,(g) =iV0,(g) + A7V„,(g); L„, _i/2 = i° , _i/2 + AL,, _i/2 . (59) 

From use of the ground-state distribution and number values given in Eqs. I|A12|I - I|A14() of Appendix A, we find 
the following operational relations, 

: N,,iq) := N„,{q) ; ly > ; : 7V,,(q) 7V,,(g) ; iy>l; : La, -1/2 := ^a, -1/2 • (60) 

Relations (j6()|l are justified by the absence of the corresponding quantum objects in the initial ground state. 

B. THE NORMAL-ORDERED PSEUDOPARTICLE MOMENTUM AND ENERGY FUNCTIONALS 

Let us introduce the momentum and energy functionals associated with the Hilbert subspace spanned by excited 
states described by small deviations (|58|l . The bare-momentum distribution deviations of these states involve a small 
density of pseudoparticles, —1/2 Yang holons, and —1/2 HL spinous. In the case of momentum, such a functional 
is obtained directly from use of the operator expression (|34|l . Introduction of the corresponding normal-ordered 
expression and replacement of the bare-momentum number operators and holon number operators by their eigenvalues 
leads to the following momentum functional, 

ito 00 +qBu 00 +<?c. ,^ 

AP= J2 ^Nco{q)q + Yl E ^NsAq)q + J2 E AA^,,(g) [(1 + i.)- - g] + - AL,, _i/2 • (61) 

The momentum functional (|61() is linear in the bare-momentum distribution function deviations for all excited states 
independently of whether the values of these deviations are small or large. Thus, this expression is valid even for large 
values of these deviations. 

As for the Fermi- liquid quasiparticles |42l |43| , while the momentum functional H61|) is linear in the bare- momentum 
distribution function deviations, the corresponding energy functional includes linear, quadratic, and higher-order 
terms in these deviations. Such an energy functional is derived by solution of the integral equations 1)19 11 - 1)21(1 for 
bare-momentum distribution functions of the general form (|59f) . This leads to deviation expansion expressions for the 
rapidity-momentum and rapidity functionals in terms of the set of bare-momentum distribution function deviations 
ANco{q), {ANcu{q)}, and {ANsu{q)} where — 1,2,.... According to Eqs. the energy spectrum of the 

excited states can be expressed in terms of these functionals. Use in the general energy expressions (I31|l - (|33|) of the 
systematic expansion in the pseudoparticle bare-momentum distribution deviations and holon and spinon number 
deviations of the above functionals leads to a general finite-energy Landau-liquid energy functional of the following 
form. 



00 

AE^Y^ AEi 

i=l 



(62) 
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The terms of order i larger than one describe the residual interactions of the pseudoparticles. In contrast, the 
energy functional H62I) is linear in the —1/2 Yang holon and —1/2 HL spinon number deviations for all excited states 
independently of whether the values of these deviations are small or large. Such behavior follows from the non- 
interacting character of these quantum objects |3n| |. In the case of excited states whose bare-momentum distribution 
function deviations involve a small but finite density of pseudoparticles, the physics is described by the energy terms 
on the right-hand side of Eq. (|62() of all scattering orders i = 1,2, .... However, in the limit when the density of these 
pseudoparticles becomes vanishing, only the first two terms become relevant. According to the results obtained in 
Refs. |2Es.30j, the first and second-order terms on the right-hand side of Eq. (|(j2|l arc of the following general form, 

ASi=c^O+ dqecoiq)AN,o{q)+ ^ e,i (g) A7V,i (g) + ^ ^ dq e^Jq) AN^^^iq) , (63) 

q=-q"g q=~q"i "=c, s ^ q=_gO_^ 

and 



oo +9L oo +9°',,.' 

^^^^lH E EE E E l^Uu..'u'{q,q')ANUq)^N^'..'{q'), (64) 

where ANc{q), ANc^^q), and AN^uiq) are the pseudoparticle bare- momentum distribution function deviations given 
in Eq. (|58fl and ujq is an energy functional linear in the —1/2 holon, —1/2 spinon, and si pseudoparticle number 
deviations AM^,. _i/2 AAfj, _i/2, and AA^^i, respectively. (The latter functional is given below in Eq. H1U1|) of Sec. V, 
with the si pseudofermion number deviation equaling the corresponding si pseudoparticle number deviation.) The 
coefBcients of the i — 1 linear terms are the af pseudoparticle energy bands ecQ{q), esi{q), and ^ai^iq) studied in Refs. 
[29! l3[ll|ul |. The coefficients of the i = 2 quadratic terms are the pseudoparticle / functions faiy,a'v'iq,q') studied in 

In contrast to what the energy terms (|63l) - (|64|l may suggest, we emphasize that the energy (I62|) is not an expansion 
in 1/L. The small parameters in such an expansion are the pseudoparticle bare-momentum distribution function 
deviations, as in the case of the quasiparticle Fermi-liquid energy functional [4^ l43l] . Indeed, when these function 
deviations involve a small but finite density of pseudoparticles, the energy terms of order i on the right-hand side of 
Eq. H62() are not of order [l/i]\ Instead, in that case all terms of the energy expansion H62() are of the same order 
[1/L]^^ — L. As discussed in the ensuing section, few-electron excitations are associated with excited states generated 
from the ground state by a finite number of pseudoparticle, holon, and spinon processes. The deviations of these 
states involve a vanishing density of pseudoparticles, holons, and spinous. However, in order to achieve the correct 
microscopic description of the few-electron spectral properties, one must consider the limit where a vanishing density 
of pseudoparticles is approached rather than considering that such a density is zero from the very beginning. As the 
density of pseudoparticles involved in the bare-momentum distribution deviations decreases and approaches zero, the 
energy terms on the right-hand side of Eq. Ht)2|) of large scattering order i, also vanish and do not not contribute as 
they become of order [1/i]'. 

The / functions on the right-hand side of Eq. (|64|l . are associated with the two-pseudoparticle residual interactions 
and have the same role as those of Fermi- liquid theory [i^lisj . Indeed, for small values of the energy lo and electronic 
densities n and spin densities m such that < n < 1/a and < m < n, respectively, the few-electron spectral 
properties are controlled by the residual two-pseudoparticle interactions described by the i — 2 terms (|64|l of the 
energy functional Ht)2|) . In the corresponding low-energy Hilbert subspace, only the cO and si pseudoparticle branches 
have finite occupancies. In this case, as the limit of vanishing density of pseudoparticles is approached, the general 
energy functional Ht)2|l acquires the form of the energy spectrum of a two-component cO and si conformal field theory 
|2l|. For the relation of the c = c, and s = s, 1 pseudoparticle description to two-component conformal- field theory 
see Refs. The weight distribution of the few-electron spectral functions can in this case be studied by 

conformal field theory [2^ . 

We thus conclude that the pseudoparticles have residual interactions associated with the terms of the general energy 
functional (|62|) of order i > I. The pseudoparticle residual interactions control the low-energy few-electron spectral 
properties of the quantum liquid, as confirmed by the studies of Refs. |2^ . In this paper, we introduce an 
alternative representation in terms of non-interactiong pseudofermions. For such pseudofermions the bare momentum 
is replaced by a momentum functional whose coefficients are two-pseudofermion phase shifts. In this case, the relevant 
energy spectrum is non-interacting and thus has only linear terms in the momentum distribution function deviations. 
For the pseudofermion description the few-electron spectral properties are controlled by such two-pscudofcrmion phase 
shifts ,34iJ. Such shifts are associated with first-order momentum distribution function deviations only, whereas the 
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pseudoparticle residual interactions involve scattering orders i associated with terms of the energy functional l|62|) 
such that i > 1. 



C. THE CONCEPT OF J-CPHS ENSEMBLE SUBSPACE 



Let us introduce the quantum number i — sgn((7)l ~ ±1 which refers to the number of right (i = +1) and left 
(l ~ —1) pseudoparticle movers. The number A^qj/, t of pseudoparticles of l character is a good quantum number. 
Wc thus introduce the av pseudoparticle current number Jai, such that, 



Each CPHS ensemble subspace contains in general several subspaces with different values for the sets of current 
numbers Jco and {JaiA such that a = c,s and u = 1,2, .... We call these subspaces J-CPHS ensemble subspaces. For 
a given J-CPHS ensemble subspace one can introduce the J-CPHS ground state. This energy eigenstate has compact 
bare-momentum occupancies for the av pseudoparticle bands. Moreover, a partial J-CPHS ground state is a state 
with such a type of occupancy configuration for the av pseudoparticle branches such that av ^ cO and av ^ si 
only. These concepts are associated with the right (t = +1) and left (i = — 1) Fermi points, which separate the av 
pseudoparticle occupied from the unoccupied regions of bare-momentum space. Ignoring terms of order 1/L, these 
Fermi bare-momentum values read. 



27r 

qFaiy,L^ l-^Nai,. ,. (66) 

The ground state is a particular case of J-CPHS ground state such that the Fermi points of Eq. read. 



'lFau,i — '-'fpav 5 (67) 

where 



9fc0 = 2fcF ; qpsi = kpi ; Qfc^ = 0; qp^^ = , ly > Sa 



(68) 



In most situations one can disregard the 1/L corrections and use the bare-momentum Fermi values provided in Eqs. 

The pseudofermion description introduced in the ensuing section, corresponds to J-CPHS ensemble subspaces 
spanned by energy eigenstates differing from the ground-state by the occupancy configuration of a finite number of 
av pseudoparticles, —1/2 Yang holons, and —1/2 HL spinous. For such subspaces the pseudoparticle Fermi band 
momenta given in Eq. H66|) can be written as follows, 



qFaiy, I. — L q%au + ^qPau, l 

where q^^j^ is the av ground-state Fermi bare momentum whose values are given in Eq. (|68|l and thus. 



(69) 



qpcu, L = Aqpcu, I, v>Q; qpsv, u = ^qpsv, l, v > 1 . 



(70) 



Here 



27r ^ 2tt 



iAJai. 



(71) 
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V. THE PSEUDOFERMION DESCRIPTION 



In this section we introduce the pseudofermion operational description. This includes introduction of the pseudopar- 
ticle - pseudofermion unitary transformation and of the Hilbert subspace associated with that transformation. In such 
a subspace, the rapidity functionals have a particular form. Moreover, we obtain the anticommutation relations of 
the pseudofermions. 



A. THE FUNCTIONAL CHARACTER OF THE PSEUDOFERMION MOMENTUM 

The av pseudofermion carries momentum qj given by, 

q, = q{q,) ^ q, + ^ ^I"-^ + ; J = 1, 2, TV*, . (72) 

Here QaviQj)/ L is the momentum functional, 

^^ = fE E E''J'-."V'fe,<Z,OAA^«'.'feO. (73) 

The functions ^af, aV (q, <?') on the right-hand side of this equation are the two-pseudofermion phase shifts expressed 
in terms of the bare-momentum q. These phase shifts are given by, 



where 5>q,^_ q/^/ (r, r') are the two-pseudofermion phase shifts expressed in terms of the variable r. The latter phase 
shifts are the unique solutions of the integral equations (|B2|) - (|B15|) of Appendix B. The value of the momentum 
functional given in Eq. H73() is controlled by shake-up two-pseudofermion phase-shift processes resulting from the 
changes in the quantum-object occupancy configurations described by the deviations ANa'u' (Qj' ) on the right-hand 
side of that equation. 

The av pseudofermion is related to the corresponding av pseudoparticle by a mere unitary transformation involving 
the discrete bare- momentum values qj and such that. 



Qj , (75) 

where qj is the discrete momentum defined in Eq. H72|) . The number of av pseudoparticles, iVo,^, equals that of av 
pseudofermions, which we denote by Afai^. Moreover, we introduce the ai^ pseudofermion momentum distribution 
function Nav{<ij) such that Afa,A<lj) — ^ai'{<lj{<ij)), where qj — qj{qj) is the inverse function of H72I) . 

Note that the momentum functional ||72J), is such that Qa^u{<l) / L = for the initial ground state. Indeed, the 
pseudofermion description refers to the ground-state normal ordered ID Hubbard model. Thus, there is a specific 
av pseudofermion description for each initial ground state. For the latter state the bare momentum qj = [27r/L]/"'' 
equals the pseudofermion momentum qj — qj + Qauiqj)/ L. This justifies the designation bare momentum for qj. 
Thus, in the case of the ground state the av pseudoparticles are identical to the av pseudofermions. It follows that 
the ground state is invariant under the pseudoparticle - pseudofermion unitary transformation and plays the role of 
the vacuum of the pseudofermion theory. 

As for the case of the pseudoparticle representation, the pseudofermion description cor resp onds to large values 
of the Hubbard chain length L and is thus compatible with Takahashi's string hypothesis [2^|31|- However, while 
the pseudoparticle representation corresponds to the whole Hilbert space, the pseudofermion description refers to a 
Hilbert subspace spanned by the initial ground state and the following types of states: 

(A) - Excited states generated from the initial ground state by a finite number of pseudofermion processes involving 
creation or annihilation of cO and si pseudofermions, creation of pseudofermions belonging to other av ^ cO and 
av 7^ si branches, and creation of —1/2 Yang holons and —1/2 HL spinous. In the thermodynamic limit, the 
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momentum and bare-momentum distribution function deviations associated with such excited states obey the following 
relations, 

ALa,_i/2/i La,-i/2/L ^ 0; a = c,s. (76) 



(B) - Excited states generated from the initial ground state by a finite number of cO and si pseudofcrmion particle- 
hole processes. 



Throughout this paper, we call such states, excited states A and excited states B, respectively. We emphasize that 
in general the excited states involve both. However, for few-electron excitations which do not change the electronic 
numbers the low-energy states are of type B only. Moreover, for simplicity we are considering densities in the ranges 
< n < 1/a and < m < n. Our analysis also holds for other values of the densities, yet for the half-filling n = 1/a 
or zero-magnetization m = phases the excitation subspace is more reduced. 

The unitary transformation associated with Eq. H75() , is well defined for the Hilbert subspace spanned by the initial 
ground state and the excited states A and B. Fortunately, in the thermodynamic limit application onto the ground 
state of a few-electron operator generates excited states A and B only. By few- electron operators, we mean here 
operators which can be written as a product of Af electron creation and/or annihilation operators and J\f/Na — > 
as Na oo. The pseudofermion functional theory introduced in this paper, can be applied to the study of the 
few-electron spectral weight distributions for all values of energy [l|, [l^, ■ We find below that the pseudofermions 
have a non-interacting character. This allows the evaluation of few-electron spectral functions for all values of energy 
[let l34l|. The ai^ pseudoparticle residual interactions are cancelled by the momentum transfer term Qa^vio) / L of 
Eq. (|73|l . The information recorded in the pseudoparticle interactions is transferred over to the momentum two- 
pseudofermion phase shifts of the momentum term Qaviq) / L. The pseudofermion theory is of first order in 1/ L both 
for the energy and momentum spectra. In contrast, the methods that use the pseudoparticle residual interactions lead 
to expressions for these spectral functions which are valid for small values of the energy lo only 0, (23. l29l | . Therefore, 
the pseudofermion description introduced below corresponds to a breakthrough in what the study of the finite-energy 
spectral properties of the ID Hubbard model is concerned. 

That only excited states A and B have finite overlap with few-electron excitations, is confirmed by direct evaluation 
of matrix elements between few-electron excitations and excited states (34j . For instance, if Qau{Q)/L — for the 
excited states B, the contribution to the one-electron spectral function would come from one particle- hole pseud- 
ofermion processes only, and the spectral functions would have the familiar (5-function structure. This is not true for 
our case where Qav{<l)/ L ^ for these states. In this case, there are contributions from many particle-hole processes 
in the cO and si pseudofermion bands as well. The largest weight comes from the one particle-hole pseudofermion 
processes, and increasing the number of particle-hole processes, the additional weight decreases rapidly. Most of the 
weight associated with the excited states B is generated by one, two, and three particle-hole pseudofemion processes 
[3^ . In the general case of few-electron spectral functions, the contribution from excited states B involving an infinite 
number of particle-hole processes in the cO and si pseudofermion bands vanishes in the thermodynamic limit. The 
same holds for the excited states A. The excited states generated from the ground state by an infinite number of 
processes of the types reported in (A) and (B) have vanishing overlap with few-electron excitations as L — > 00. 

The cancellation of the av pseudoparticle residual interactions by the momentum transfer term Qa.vio) / L of Eq. 
(|73ll is related to the form of the rapidity functionals Ka^^q) and rapidity-momentum functional k(q) in the Hilbert 
subspace spanned by the excited states A and B. Introduction of the pseudoparticle bare- momentum distribution 
functions of general form given in Eq. I|59|) in the rapidity functional integral equations (|19l) - (|21|l and their expansion 
in the small deviations (|58|l . permits explicit solution of these equations. This procedure leads to expressions for 
the rapidity functionals Kay{q) and rapidity-momentum fimctional k{q) in terms of the bare-momentum distribution 
function deviations introduced in Eq. (|58|l . Solution of the integral equations H19|) - (|21|l for distributions of the general 
form H59|l leads to first-order in the deviations to expressions for the rapidity-momentum functional and rapidity 
functionals of the following form. 



k{q)^k°[q{q)y, A„,(g) = A" , (^(q)) ; a = c, J. = 0, 1, 2, ... ; a = s, z/=l,2,.... (77) 
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Here hP^^{q) and k'^{q) are the corresponding ground state functions. On the right-hand side of Eq. (|77|) . q{q) is the 
av momentum functional given in Eq. (|72() with qj replaced by the continuum momentum q. 

It is remarkable that in the Hilbert subspace spanned by the excited states A and B the functionals h-auio) and k{q) 
equal the corresponding ground-state functions A^j,(q) and k^{q), respectively, with the bare momentum q replaced by 
the momentum functional (|72|l . This property is behind the non-interacting character of the pseudofermions. The two- 
pseudofermion phase shifts $q^, qvC?", 't') defined by Eqs. ljB2|) - (|B15|) of Appendix B which appear in the expression 
of the momentum functional H73(l through the relation (|74|l play a central role in the pseudofermion description of 
the few-electron spectral properties. Within that description the expression of many physical quantities involve such 
two-pseudofermion phase shifts. For instance, the functions defined by Eqs. (|47l) - (|49|l and the expressions of the 
pseudofermion energy bands given below involve the two-pseudofermion phase shifts. Furthermore, these phase shifts 
control the matrix elements between few-electron excitations and the excited states and thus also the few-electron 
spectral weight distributions [33 |. 

According to Eq. H33() the energy spectrum depends on the quantum object occupancy configurations through 
the rapidity and rapidity-momentum functionals. It is this functional character that is behind the pseudoparticle 
residual interactions and leads to the energy expansion given in Eqs. H62(l - (|64|l . However, by re-expression of these 
functionals in terms of the pseudoferminon momentum q, the general energy spectrum H31|) - (|33[l can be written in 
terms of pseudofermion momentum distribution functions McoiQj) and Mciy{qj) as, 

Na OO N'^ 

E - -2t£ AA,ofe) cosfc"fe) +4t^^ AA,,fe)Re{y/l - {M^.^q^) + ivU / Uy] 



U 



3 = 1 u=l j=l 

N„ 



M,^J2^iyK.-^\+J2tiaS^. (78) 



This expression is valid for the Hilbert subspace spanned by the initial ground state and the excited states A and 
B. The term t^aS" is the same as on the right-hand side of Eq. (|31l) and Afd, = N^, is the number of ciy 
pseudofermions. The form of the general energy spectrum (|78|l justifies why the shake-up effects associated with 
the phase shifts of the functional ]T6\i occur in the case of the pseudofermions in the momentum instead of in the 
energy. The dependence of the general energy spectrum (|78|l on the momentum occupancy configuration values occurs 
through the arguments of the ground-state rapidity and rapidity-momentum functions. Thus these functions play the 
role of non-interacting spectra, since they have the same form both for the initial ground state and excited states A 
and B. The shake-up effects associated with the two-pseudofermion phase shifts are thus felt by the non-interactiong 
pseudofermions as mere changes in the momentum occupancies through the twisted boundary conditions imposed 
by the elementary processes which generate the excited state from the ground state. This property is behind the 
non-interacting character of the pseudofermions. It is also behind the fact that within the pseudofermion description 
the functions qdk) and qa,y{A) defined in Eqs. H47l) - (|49|1 refer to both the ground state and the excited states. In 
contrast, within the pseudoparticle representation these functions refer to the ground state only. 

The pseudoparticle bare-momentum q, description is naturally provided by the Bethe-ansatz equations within 
Takahashi's string hypothesis [2^ l29l l31|. We recall that the pseudoparticle discrete bare-momentum values qj are 
of form given in Eq. ljA2p of Appendix A and according to Eq. (jAl|) are such that g^+i — qj — 2-k/L. The single 
discrete bare-momentum values qj are integer multiples of 2'k/L or of ir/L jsj and bare- momentum contributions 
of order [l/L]^ such that j > I are outside the vahdity of the pseudoparticle description. These bare-momentum 
contributions have no physical meaning and must be considered as zero. Importantly, the same is required for the 
pseudofermion momentum discrete values qj given in Eq. \I2\ . These discrete values are also at least of order of 1/L 
and contributions of order [1/i]-' such that j > 1 must be considered as equaling zero. For instance, as confirmed in 
Appendix C, the discrete momentum level separation, 



27r Qai/(gj-i-i) - 27r 

~ ^ "/7 ^ 1 ]7 ' ^ '^1 

is such that the second term on the right-hand side of Eq. 1)79(1 is of order [1/L]^, where Qau{q)l L is the momentum 
functional given in Eq. I(73() . Thus, to first order in 1/L one finds that qj+i — qj = 2'k/L, as for the corresponding 
discrete bare-momentum level separation given in Eq. I|A1|I of Appendix A. However, this does not imply that to first 
order in 1/L the pseudofermion momentum equals the bare- momentum. Indeed, we emphasize that the momenta 
Qau{qj)/ L on the right-hand side of Eq. 1(75(1 are of order 1/L and play a central role in the control of the few- 
electron spectral weight distribution by the non-perturbative many-electron shake-up effects We note that the 
level separation qj+i — qj = 2'k/L is valid locally in the discrete momentum space. By that we mean the following: If 
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in the present thermodynamic limit two momentum values qj and qji differ by a small yet finite momentum difference 
= qj — (jji, then in general Aq ^ ^[j — j']. In contrast, for the corresponding bare- momentum values it holds that 
Aq — ^[j — j']- Therefore, for small but non-vanishing momentum separation the difference [Qai^{qj) — Qawiqj')]/ L 
is not anymore of order [1/ LY and has physical significance. 



B. PSEUDOFERMION OPERATOR ALGEBRA 



The elementary creation and annihilation operators of the av pseudofermions can be expressed in terms of the 
corresponding operators of the av pseudoparticles as follows, 

f\j , av ~ '^av , au ^av ] fqj , av — ^au ^<3j , au ^av ■ (^0) 

Here Vaii is a unitary operator that we call the av pseudoparticle - av pseudofermion unitary operator. This operator 
acts onto the Hilbert subspace spanned by the initial ground state and excited states A and B. It shifts the bare- 
momentum value qj by Qai, {qj ) / L and thus is given by, 

VL = exp{^ ^g,+Q_fe)/L,a.^9.,"-} 5 = ^Xpj^ ^9, (9, )/L. a,. , } " (^1) 

The momentum distribution function ^>fau{qj) is the eigenvalue of the operator, 



J^av{qj) = fq^^au fq,.av ] J^av{qj) = Nau{q]{qj)) ■ (82) 

Here Nav{qj) stands for the av pseudoparticle bare-momentum distribution function and qj{qj) is the inverse of the 
function H75|l . The function Na^{qj) is the eigenvalue of the corresponding bare- momentum number operator H18|l . 
Keeping only the physical contributions that correspond to terms up to first order in 1/L, the function qj{qj) is given 
by, 



<l:=lA^,) = 'ij^^^-^,-Y E E E C'.'fe-,9.')AAC.,.,feO- (83) 

a'— c,s v'—h^i ^ j'—^ 

Here, AA/'qV'(9j') is the deviation of the function (|S^ relative to its ground-state value. Since the functional l(75|l 

vanishes for the ground state, note that qj = qj for that state. Moreover, on the right-hand side of Eq. H83|l 
^ai' 9') is the two-pseudfermion phase shift. It is defined as, 



= *.....(*),*-')) = t...v. , (M) 

where the function q = q{q) is the continuum version of the function H83|l . ^av.a'v'{q,q') is given in Eq. H74II . and 
^afja'v'if, r') is the two-pseudofermion phase shift expressed in the variable r defined by the integral equations 
(|B2l-(IHT5j of Appendix B. 

Often in the expressions given in previous sections we replaced the pseudoparticle bare-momentum summations by 
integrals and the corresponding discrete bare-momentum values qj by a continuum bare-momentum variable q. Since 
according to Eq. IjAip of Appendix A, the difference g^+i — qj = 27r/L is constant for all values of j, the use of such 
a continuum representation involves the replacement. 



In the Hilbert subspace spanned by the excited states A and B, the rapidity functional Aa^{q) and rapidity- 
momentum functional k{q) equal the corresponding ground-state rapidity function Aj^^(g) and rapidity- momentum 
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function k^{q), respectively, with the bare- momentum q replaced by the momentum q. It follows that in such a 
subspace the limiting values of the continuum momentum q are given by the ground-state limiting values igj^^ given 
in Eqs. H39|l and (|4()|l . (This is confirmed in Sec. Vll.) Thus, to replace the discrete momentum values by a continuum 
momentum variable q, one must replace the summations by the following integrals. 



2tt 



dq 



dq{q) 
dq 



(86) 



We then introduce the momentum distribution function. 



X,,(g) = ^A/;,.(g). (87) 
dq 

The function q = q{q) on the right-hand side of Eq. H86|) is the inverse of the function given in Eq. (|72|l and is given 
in Eq. H83(l with qj replaced by q. Its derivative is given by, 



-^ = l-Y. Ey dq AAA^vK-z). (88) 

The second term on the right-hand side of this equation is of first order in the momentum distribution function 
deviations. For the non-interacting pseudofermion description only momentum and energy contributions up to first 
order in these deviations are physical. As a result, in the case of momentum distribution function deviations ^Mavisi) 
one can consider that. 



Ayv;,(g) = AAA„,(g), (89) 



where in contrast to the case of Eq. H87|l we used dq{q)/dq — 1. 
The av pseudoparticle number operator. 



j=i q=-qa 



"q., av "9: av 



2lT 



dqh\ 



q, av "q, av , 



(90) 



is invariant under the pseudoparticle - pseudofermion transformation. It equals the number of av pseudofermions 
operator, 



+qc. 



J^av ^ ^ fqj,avfqj,av ^ ^ fq,avfq-av 



27r 



,- dq{q) I 
dq ""^ 



/q, av ■ 



(91) 



The cv pseudoparticle charge and Sc and Cc values considered in Sec. Ill are also invariant under the same transforma- 
tion. The same occurs for the sv pseudoparticle Sg and Ug values. The pseudoparticle - pseudofermion transformation 
also leaves invariant the ±1/2 Yang holons and ±1/2 HL spinous. The ±1/2 holon (and ±1/2 spinon) composite 
character of the cv pseudoparticles (and sv pseudoparticles) also remains invariant under that transformation. It 
follows that the cv pseudofermions (and sv pseudofermions) are Sc = (and Ss = 0) composite objects of an equal 
number v = \,2, ... of —1/2 holons and +1/2 holons (and —1/2 spinous and +1/2 spinous). Consistently, the ±1/2 
holon (a = c) and ±1/2 spinon (a s) number operators Af^, -1-2/2 given in Eq. (|29|l can be rewritten in terms of 
pseudofermion operators as follows. 



±1/2 



±1/2 



E 

v=l 



+9L 

E 

q=-ql 



vNav{q) : 



(92) 
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where the pseudfermion momentum distribution operators A/'ay (5) are given in Eq. (|82l) . Like in Eq. 129|) . the operator 
-£a,±i/2 on the right-hand side of Eq. H92I) is the ±1/2 Yang holon (a = c) and ±1/2 HL spinon {a — s) number 
operator given in Eq. (|30|l . AU results reported in Sec. Ill concerning pseudoparticle charge and spin transport are 
also valid for the corresponding pseudofermions. For instance, for finite values of U/t the transport of charge (and 
spin) is associated with the cO pseudofermion and cv pseudofermion quantum charge fluids (and sv pseudofermion 
quantum spin fluids). 

It is found in Ref. and discussed in Appendix A that the bare-momentum q is the conjugate of the spatial 
coordinate Xj — aav j associated with the effective av lattice, where j = 1,2, ...,N*j^. As in the case of the charge 
(or spin) carried by the pseudoparticles and of their composite character in terms of chargeaons and antichargeons, 
±1/2 holons, or ±1/2 spinous, also the effective av lattice remains invariant under the av pseudoparticle - av 
pseudofermion unitary transformation. Indeed, the momentum functional Qau{q)/L which controls the pseudoparticle 
- pseudofermion transformation H75|) does not affect the underlying effective av lattice. That momentum functional 
just imposes a twisted boundary condition such that each av pseudofermion hopping from site N*^^ — 1 to site of 
such an effective lattice will acquire a phase e^'^a,.!'?)^ From combination of that analysis with the expression for the 
general momentum functional Qau{q)/L given in Eq. I|73|) . we find that the phase shift <I>^^ av(9'9') such that 
creation into the system of a a'v' pseudofermion of momentum q' imposes a twisted boundary condition on the wave 
function of the av pseudofermion branch. This condition is such that when the av pseudofermion of momentum q 

hops from site iV*^^ — 1 to site of the effective av lattice, it will acquire a phase e*^^*"" "'"' The studies of 
Ref. |3| reveal that within the pseudofermion description the few-electron spectral properties are fully controlled by 
such two-pseudofermion twisted boundary conditions associated with the phase shifts 

As for the case of the pseudoparticles, it is useful to introduce the local av pseudofermion creation operator f^. 

and annihilation operator fxj,av These operators are related to the operators f^q^av ^-^d fq,av, respectively, obtained 
from the corresponding pseudoparticle operators through the relations given in Eq. H8U|I. as follows. 



n = — 



,iqxj p 



fq. av /— 



J Xj , av 7 



(93) 



where the summations refer to the sites of the effective av lattice. The local av pseudofermion creation (and anni- 
hilation) operator f^. (and fxj,av) creates (and annihilates) a av pseudofermion at the effective av lattice site 
of spatial coordinate Xj = a'^^j, where j = 1,2, N*^ and a^j, is the effective av lattice constant given in Eq. 
(|37|l . (Since the pseudofermion representation refers to the Hilbert subspace spanned by the initial ground state and 
the excited states A and B, except for 1/L corrections we can consider that the corresponding effective av lattice 
constants are the ground-state constants a^^.) Thus, the conjugate variable of the momentum qj of the av pseud- 
ofermion branch is the space coordinate Xj of the corresponding effective av lattice. The local av pseudoparticles and 
corresponding local av pseudofermions have the same effective av lattice. It follows that the local pseudoparticle and 
local pseudofermion site distribution configurations which describe the ground state and the excited states A and B 
are the same. These configurations are expressed in terms of rotated-electron site distribution configurations in Ref. 

While the local av pseudoparticles and corresponding local av pseudofermions live in the same effective av lattice, 
the values of the set of discrete bare-momentum values {qj} and momentum values {qj} such that j = 1,2, ...,N*^ 
are different and related by Eq. (|83|l . There is an one-to-one relation between these two sets of discrete values, which 
keep the same order because there is no level crossing. This property follows from the discrete bare-momentum and 
momentum separation given in Eq. (|A1() of Appendix A and Eq. (|79f) . respectively. 

Finally, we consider the anticommutation relations of the pseudofermion operators. It is confirmed in Ref. [3^ 
that such relations play a major role in the evaluation of matrix elements between energy eigenstates. Let us consider 
the general situation when the momenta q and q' of the operators f^^av ^-^d fq\av, respectively, correspond to 
different J-CPHS ensemble subspaces. The anticommutator fq',a'v'} can be expressed in terms of the local- 

pseudofermion anticommutators {f^. fxj,.a'v'} associated with spatial coordinates Xj and xj' of the effective av 
and a'v' lattices, respectively, as follows, 



{fl.avi fq',a'v'} ^ j~ ^ ] 6 ^ ^^{fxj,avjfxj,,a'v'}- (94) 
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After performing the j and j' summations we find that the av pseudofermion operators obey the following algebra, 

t 1 f /. fn M n' , sin([g„.(g)~g'„.(q')]/2) 

^' L sm([(7 — q'\ a/2) 

and the anticommutators between two av pseudofermion creation or annihilation operators vanish. On the right-hand 
side of Eq. H95|l the functionals Qav{q) and Q' civil') whose general expression is given in Eq. (|73|l refer to the J-CPHS 
ensemble subspaces where the momenta q and q' refer to, respectively. Note that when s\n{[Qav{q) — Q' auil)]/'^) 
vanishes the anticommutation relation (|95|l is the usual one, {/^^q,^, fq',a'u'} — Sa,a' Si^.u' ^q.q'- In contrast, for finite 
values of that quantity there is an overall two-pseudofermion phase shift which arises from a non-perturbative shake-up 
effect associated with the functional character of the momentum Qauiq)/L given in Eq. I|73|l . 

A case of particular importance is when one of the two J-CPHS ensemble subspaces associated with the momenta q 
and q' is that of the initial ground state. In such a case Q' avW) — ^r the ground-state J-CPHS ensemble subspace 
and thus the anticommutation relation H95|l simplifies to. 



VJ q, av 5 



-i{q- 



^q')a/2 ^iQc 



-(■z)/2 



3in(Qa.(g)/2) 
sin([g-g'] a/2) 



(96) 



VI. THE PSEUDOFERMION ENERGY AND MOMENTUM SPECTRA AND THE WAVE-FUNCTION 
FACTORIZATION OF THE NORMAL-ORDERED ID HUBBARD MODEL 

In this section, we find that the description of the quantum problem in terms of the non-interacting pseudofermions 
leads in the thermodynamic limit to a wave-function factorization for the excited states A and B contained in the 
excitations generated by application onto the ground state of few-electron operators. Such a factorization refers to 
the ground-state normal-ordered ID Hubbard model. 



A. THE PSEUDOFERMION ENERGY AND MOMENTUM SPECTRA 



Provided that all scattering orders are considered, the pseudoparticle energy functional Ht)2|l describes excitations 
which involve a small, but finite density of pseudoparticles. However, in the thermodynamic limit, the few-electron 
excitations are contained in a smaller Hilbert subspace. This subspace is spanned by the excited states A and B. 
These states are generated from the ground state by processes which involve a finite number of pseudofermions (and 
pseudoparticles). In the thermodynamic limit, this corresponds to a vanishing density of these objects rather than to 
a small but finite density. In the Hilbert subspace spanned by the ground state and the excited states A and B the 
energy spectrum is of the form given in Eq. (|78|l . This spectrum can be expressed in terms of continuum momentum 
integrals with the result. 



E = -2t— / dgAr,o(g) cosfc"(g)+4t— V / dqMc.isi) Rc | ^1 " {^"iM + ii'Vl^tf \ 



U 



Na 



where Mco{q) and Afci^{q) are the pseudofermion momentum distribution functions of Eq. I|87|l . 

In contrast to the pseudoparticle energy functional (|62|l -l|64 () in terms of the pseudoparticle bare-momentum distri- 
bution function deviations, the energy spectrum derived from expressions (f75|) and only includes first-order terms 
in the pseudofermion momentum distribution fimction deviations. In the Hilbert subspace that the energy spectrum 
(|78|l - (|97|l refers to, only such first-order pseudofermion momentum distribution function deviation contributions have 
physical meaning. The information recorded in the pseudoparticle interactions is transferred over to the momentum 
two-pseudofermion phase shifts of the momentum term Qai,{q)/ L. The pseudofermion discrete momentum values are 
of 1 /L order and store this information. 
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According to Eq. (|73|) , the processes which generate the excited states A and B from the initial ground state lead to 
a collective momentum shift of all the cO pseudofermions and si pseudofermions of the initial-state Fermi sea. That 
momentum shift reads, 

+ = 0, ±7r; a = cO, si, (98) 

where the bare-momentum shift t:^^/ L is given by, 

oo oo 

T^co = 0; ^l^Na,^ even; 7r°o = ±7r ; ^ ^ANa^ odd; 

Q— c, s v—l OL—c. s y—1 

7r°i = ; AiVco + AiV^i even ; 7r°i = ±7r ; AiV^o + AA^si odd . (99) 

The excited states A and B are generated from the ground state by two virtual excitations: (1) a collective momentum 
shift of the Fermi sea cO pseudofermions and si pseudofermions; (2) for the states A, a finite number of cO pseud- 
ofermion and si pseudofcrmion creation or annihilation processes and of av ^ cO and av ^ si pseudofermion, —1/2 
Yang holon, and —1/2 HL spinon creation processes; for the states B, a finite number of cO pseudofermion and si 
pseudofermion particle-hole processes. Consideration of both the energy associated with the virtual excitations (1) 
and (2), gives the energy spectrum of few-electron excitations associated with the energy (|78|l - (|97|) . Such an energy 
spectrum corresponds to the ground-state normal-ordered ID Hubbard model and is additive in the —1/2 holon, —1/2 
spinous, and av pseudofermion energies. It reads. 

No. N*^ oo N*,^ 

AE = ^0 + ^ AKo {qj ) eco {Qj ) -f £ AA^.i (g, ) e,, i (g, ) -f ^ ^ g AA^., (g^ ) 6° , (g) 

j=l j = l a=c, s i/=l+S^ ,. j=l 



= t^O + 

L 



— / dqAAfrMiq)ecoiq) + 7rz2 dg AAA,i(g) e,, i(g) 



' 9 E E / dqAM^,iq)elM, (100) 

a = C., S iJ=l + 5a,, a 

where the energy parameter cjq is given by, 

ujQ = 2^l AA4, _i/2 + 2/10 H [AM,, _i/2 - AA^.i] . (101) 

Here AMq, _i/2 are the deviations in the numbers of —1/2 holons (a = c) and of —1/2 spinous (a = s) and 
AMsi = ANsi is the deviation in the number of si pseudofermions. The limiting-momentum values gJ?o, ^sIj ^^'^ q'av 
appearing in the integrations of Eqs. (|97|) and l|100|l are given in Eqs. H39|l and (|40|l . Moreover, on the right-hand 
side of Eq. H100|l the functions eco(g), es,i(g), and ea^(g) are the pseudofermion energy bands. These energy bands 
equal the corresponding pscudoparticle energy bands appearing on the right-hand side of Eq. (|63|1 . provided that the 
bare momentum g is replaced by the momentum g. The latter bands are plotted in Figs. 6 to 9 of Ref. as a 
function of the bare momentum for zero spin density and several values of electronic density and on-site Coulombian 
repulsion. 

The shape of the pseudofermion energy bands is controlled by the two-pseudofermion phase shifts and reads, 

U ~ 
eco{q) = -^-'2tcosk°{q)+2t dk^^o.co ik,k°{q)) sink + fi - fioH , (102) 

2 J-Q 



e°,(g-) = -^C/ + 4tRc{yi - [AO,(g) - zz.^]2} + 2t j''\k ^.o^cA^, ^cM) sink, (103) 
eliq) = 2t / &$,o... (fc, A°,(g)) sinfc, (104) 

J-Q 
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and 

e.i(g) = 4(9) + 2Moi^. (105) 

The two-pseudofermion phase shifts on the right-hand side of these equations are defined in Eqs. (|50|I - H51|) . The 
parameter Q is given in Eq. H45() and the functions k'^{q) and A°^(q) are the inverse of the functions defined in Eqs. 
(|47|l and l|48(l . respectively. The zero-energy levels of the energy dispersions H1Q2(I - (|1()5|I are such that, 

eco(±2fcF) = e.i(±fcFi) = eU±[^/a - 2kF]) = e^Mkr^ ~ kpi]) = . (106) 

The dependence on the electronic density n, spin density to, and U/t of the chemical potential /i and magnetic field 
H appearing in Eqs. H1U2|I and H1U5|I can be expressed in terms of the two-pseudofermion phase shifts as, 

fi= V2tcosQ-2tl dk^cO,co{k,Q) sink - t dk ^^o. si {k,B) sink , (107) 

2 J-Q ' J-Q 

and 

t ~ 

H= / dfc$cO si (fc,-B) sin/c, (108) 

Mo J-Q 

respectively. 

We emphasize that the virtual-excitation (1) energy spectrum vanishes and thus the general deviation- linear energy 
spectrum (|100|) amounts to the contributions from the excitation (2). By construction of the pseudofermion description 
subspace, the virtual excitation (2) of the states A and B, involves changes in the occupancy configurations of a finite 
number of quantum objects. However, the excitation (1) is a collective momentum shift of all the ground-state cO 
pseudofermions and si pseudofermions. In the thermodynamic limit, the value of the corresponding ground-state 
numbers A/'co = N and Msi = -/Vj. approaches infinity. Thus, the self-consistency of the non-interacting pseudofermion 
theory implies that the energy of the excitation (1) vanishes in that limit, so that the total energy is additive in the 
corresponding pseudofermion energies. Let us confirm that as L — > oo this holds true. The cO and cl pseudofermion 
momentum distribution function deviations of Eq. H100|) can be written as. 



AAA„,(g,) - AAA«(g,) + AAA£)(g,) ; ai^ = cO, si , (109) 

where {<lj) and AAf^J (qj) are associated with the collective pseudofermion Fermi sea excitation (1) and exci- 

tation (2), respectively. Within the first order 1/L pseudofermion description the former deviation reads, 

AAf(]){q) = <^(,+ [<ii±^£^) _<,(,) 

^ [QU<l) +-L]dM^ iq) ^ [^L + Qo.(Bgn(g-)gLJ] ^(^0 _ I -|) . ^ ^ ^1 . (110) 

L oq Na 

Here, the momentum Fermi value q^^^, of the cO and si branches is given in Eq. H68|) and tt^j^/L is the corresponding 
bare- momentum shift of Eq. (|99|l . To first order in 1/L, use of Eq. (|110l) in the energy spectrum (|10U|I leads to the 
following excitation (1) energy spectrum, 

^^co!.i= E E AAAW(g-)6„.(q-) 

ai^— cO,sl qj— — q^^ 

In order to obtain the result Hlll() we have used the symmetry e^viq) — ^av{—q) and Eq. (|106|l such that £ai^{qpav) = 
for av = cO,sl. We recall that the occupancies of other ai^ ^ cO, si pseudofermion branches vanish for the ground 
state. Since the excitation (1) involves all ground-state cO and si pseudofermions, the evaluation of the 1/L higher- 
order contributions to the corresponding energy spectrum involves non-linear deviation terms. We find that such 1/L 



Y: ^ ^"-^^J^^'^-^^ e..(gLJ-0. (Ill) 

— ±1 Ql^ — CU,Sl 
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higher-order energy contributions also vanish as L ^ oo. Thus, the self-consistency condition that the excitation (1) 
energy spectrum vanishes as the system length L approaches infinity is fulfilled. 

The energy (1101(1 on the right-hand side of Eq. ((10()|l controls the finite-energy physics. The remaining energy 
terms correspond to gapless contributions provided that the involved pseudofermions correspond to the momentum 
values of the energy-band arguments of Eq. H1U6|) . For most excited states, the latter terms also lead to finite-energy 
contributions. The typical value of such energy contributions is of the order of the pseudofermion energy dispersion 
band-width per pseudofermion involved in the excited states. We note that the energy spectrum ((1()0|I of the excited 
states A and B, can have any finite value associated with the regions of the (fc, oj) plane where the few-electron 

spectral functions have finite spectral weig ht mill. 

Also, the excitation momentum spectrum can be written in terms of the pseudofermion momentum distribution 
function deviations. It is given by. 



j=l u=l j=l i/=l J = l 



-AA4,_i/2 + — / dqAN'coiq)q 



+ / dqAAfUq)q+—y2 / dq AJ\f,,{q) [- - q,] . (112) 

The important point is that the large-L pseudofermion operational description introduced in this paper, associated 
with the momentum distribution function deviation first-order energy spectrum (|100() and momentum spectrum (|112|l , 
contains full information about the few-electron spectral properties. Besides the expression in terms of pseudofermion 
operators of the diagonal (in the basis of the energy eigenstates) operators associated with these spectra, such a 
description is suitable for the evaluation of few-electron spectral functions for finite values of energy |34| . 

When acting in the Hilbert subspace spanned by the initial ground state and excited states A and B, the ground- 
state normal-ordered ID Hubbard model and momentum operator can be written in terms of pseudofermion, —1/2 
Yang holon, and —1/2 HL spinon operators as follows, 



H 



E E 



/ , ^avisij) ■ fq^.av 



■ h 



^Loi, -1/2 -^a, -1/2 : 



(113) 



a— C, S i/—5a^ s J — 1 



and 



Na 00 

j=l v=l j=l 

+ E E[(l + ^)J-9j] :/i-::/,-„c. :+-ic,-i/2, (114) 

u=l j=l °' ° 

respectively, where -/V*q = Na and the operator L^_i/2 is given in Eq. H3U|I . On the right-hand side of Eq. I|113|) . the 
cO and si pseudofermion energy bands are given in Eqs. (|102|l and (|105|l . respectively, and the av ^ cO and av ^ si 
pseudofermion energy bands and —1/2 Yang holon and —1/2 HL spinon energies read [3(tI |. 

eciy(g) = eci,(g) + 2i/^ ; es,y(g) = e"i,(g) -f 2j//io^^ ; v>5a,s, (115) 

and 



clc, -1/2 = 2/i ; eLs, -1/2 ^ 2fio H , (116) 

respectively. The energy dispersions e^^{q) and e%{q) of Eq. H115|) are defined in Eqs. H1U3|) and (|104|l . respectively. 
The ground-state normal-ordered Hamiltonian (|113|) and momentum operator 1114|) correspond to the energy and 
momentum spectra given in Eqs. H1UU|) and ()112|) . respectively. 
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B. WAVE-FUNCTION FACTORIZATION OF THE NORMAL-ORDERED ID HUBBARD MODEL 



It is well known that both the ground state wave function and the wave function of the excited states of the ID 
Hubbard model can in the U /t oo limit be constructed as a product of a spin-less fermion wave function and a 
squeezed spin wave function ^16, 18, 19J. In our pseudofermion, Yang holon, and HL spinon language this factorization 
means that in such a limit the expression of the momentum and energy of these states is linear in the av pseudofermion 
momentum distribution functions and in the —1/2 Yang holon and —1/2 HL spinon numbers. 

Let us show that for finite values of U /t the energy H31|l is not linear in the av pseudofermion momentum distribution 
functions Mavisi)- Our analysis refers to the few-electron excitation subspace. This confirms that the above type of 
factorization does not occur for finite values of the on-site repulsion. For the pseudofermion Hilbert subspace spanned 
by the ground state and excited states A and B the form of the general energy spectrum lPT|l simplifies to, 



E = E° + AE: 



E'=E% + ^[^-K.o 



(117) 



where AE is the general energy excitation spectrum given in Eq. H100|l . Eq stands for the ground-state energy, Mco 
is the number of cO pseudofermions, and the term l^aS" is the same as on the right-hand side of Eq. (|31|l . From 
combination of the ground-state occupancy configurations provided in Appendix A with the general energy spectrum 
given in Eq. we find that the energy E^^ on the right-hand side of Eq. (|117|l reads, 



E'h = -2t^ rf^ATO (g) cosA:0(g-) , (118) 

where A/Jq(^) is the ground-state cO pseudofermion momentum distribution function. (We recall that q = q for the 
ground state.) From use of Eqs. H47|) . USUI) , and H84I) we arrive to the following expressions. 



k{q)^q~J^^dq A/,o(g ) ^^q.cO (g : g) ; -^r- = I - j dq U^^iq) — . (119) 

By iterative solution of the equations given in H119|) one can derive a functional representation for the rapidity- 
momentum function kP{q) in terms of the ground-state distribution function A/Jo(9)- Use of such a functional on 
the right-hand side of Eq. H118|l leads to the following non-linear energy functional in the ground-state momentum 
distribution function M%(q), 



-2t- / dqM^,iq) 



i-iy 

m 



c«9iA(?o(9i)*Lco(9i.9) X 



OO / 



1=1 i=l 



dqi+i McQ{qi+i) 



^Ho,coi1^+l,q^) 



2i 



(120) 



Here the indices i = 1,2,... of the momentum q.i label independent integration continuum variables rather than 
discrete momentum values. The form of expression (|120() confirms that for finite values ofU/t the energy H118|l is 
highly non- linear in the ground-state pseudofermion momentum distribution function M^Q{q). Therefore, the above 
type of factorization does not occur in general for the ID Hubbard model. However, one can confirm that in the 
limit U/t oo the energy functional H120|l becomes linear in J\f^Q{q). Indeed, one finds from analysis of the integral 
equations provided in Appendix B that for spin density m — the two-pseudofermion phase shift on the right-hand 
side of Eq. H12U|I is such that ^{.q co(9'9') ^ as U/t oo. Thus, in that limit the ground-state expression H12U|I 
simplifies to. 



E% = -2t^ dqN%{q)cos{q); U/t^^. (121) 

This property also holds for the excited states and is behind the full factorization of the wave functions used in the 
studies of Refs. |16..18..19j . 
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Fortunately, the evaluation of few-electron spectral functions can be achieved without the full factorization of the 
wave functions. Such a problem can be solved by use of the ground-state normal-ordered ID Hubbard model. When 
expressed in terms of pseudofcrmion operators, that normal-ordered Hamiltonian and associated momentum operator 
are given in Eqs. and H114|l . respectively. These quantum problems correspond to the energy and momentum 

spectra given in Eqs. (|100|l and (|112|l . respectively. The ground-state normal-ordered description is equivalent to 
specifying the energy eigenstates by means of the momentum distribution function deviations AJ\fa,y{q) = AJ\fa,^(q) 
instead of in terms of the corresponding full momentum distribution functions May{q). 

In the thermodynamic limit, there is a wave- function factorization for the excited states belonging to the Hilbert 
subspace of the pseudofcrmion normal-ordered Hamiltonian H113() . This factorization follows from the form of the 
expressions (jlOOf) and (|112|l for the energy and momentum, respectively. The energy spectrum defined in Eqs. UlUUfl 
and UlUlf) can be written as. 



AS- ^ A£;„, -f 2/iAL,,_l/2 + 2Moi^AL,,_l/2, (122) 

a—c, s v—Sct s 

where the contributions from each branch arc linear in the av pseudofcrmion momentum distribution function devi- 
ations and read, 



AS^^ = V AAC.^(gj) e„^(g) = — / dqAAfa,^{q) e„^(g) . (123) 

Here the pseudofcrmion energy bands are given in Eqs. H1U2|) . (I105|l . and (|115|l . The —1/2 Yang holons, —1/2 HL 
spinous, and each pseudofcrmion branch leads to a different energy term linear in the corresponding number deviation 
or momentum distribution function deviation. We emphasize that the virtual-excitation (1) energy AE^q\^ of Eq. 
(|lll|l associated with the momentum shift of Eq. H98|l which involves all cO pseudofermions and si pseudofermions of 
the ground-state Fermi sea, is also additive in the energy contributions from each of these quantum objects. Indeed, 
according to Eqs. (|110|l and (|lll|l . the energy contribution from each cO pseudofcrmion and si pseudofcrmion involved 

in such a collective excitation vanishes. As the total energy AE^q\-^ of the same excitation also vanishes, it is indeed 
additive in the energies of all involved quantum objects. Therefore, the wave function of the energy eigenstates of the 
normal-ordered Hamiltonian can be expressed in the pseudofcrmion subspace as a product of wave functions. Each 
wave function corresponds to a different pseudofcrmion branch. In excited states with finite —1/2 Yang holon and 
— 1/2 HL spinon occupancy there is also a wave function for these objects. This factorization is associated with the 
deep physical meaning of the pseudoparticle - pseudofcrmion transformation H85|l : it transfers the information recorded 
in the pseudoparticle interactions over to the pseudofcrmion momentum, providing a non-interacting character to the 
latter objects. 

In contrast, for the pseudoparticle representation the energy functional (|62|) - (|64|l includes bare- momentum distri- 
bution function deviation non-linear terms associated with the pseudoparticle residual interactions. For instance, 
the quadratic energy term H64() contains bare-momentum summations over products of deviations of the form 
fav;a'u'{q,q') A.Nau{q) ANa'v'iq')- Such / function terms are associated with the residual two-pseudoparticle in- 
teractions. The occurrence of these energy terms mixes contributions from different branches. It follows that the 
pseudoparticle energy spectrum is not additive in the av pseudoparticle branch contributions, in contrast to the 
pseudofcrmion energy spectrum given in Eqs. (|100|l and (|122|l . Thus, in this case the wave function of the excited 
states does not factorize in the form of a product of pseudoparticle wave functions. 

The number of wave functions contributing to the factorized wave function of a given energy eigenstate depends on 
the occupancy configurations of that state. Only the av pseudofcrmion branches with finite occupancy in the state 
contribute the wave function. This contribution is in the form of a av wave function factor. The same applies to the 
occupancy of —1/2 Yang holons and —1/2 HL spinous. 

VII. TRANSFORMATION LAWS UNDER THE PSEUDOPARTICLE - PSEUDOFERMION SUBSPACE 

UNITARY ROTATION 

In this section, we find the transformation laws under the pseudoparticle - pseudofcrmion Hilbert subspace unitary 
rotation of several quantum objects and quantities and discuss the physical meaning of these laws. For instance, the 
av pseudoparticle number operator H9Uf) , is invariant under such a transformation and thus equals the corresponding 
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av pseudoparticle number operator (|91|l . Also the Yang holons and HL spinons, the charge and spin carried by the 
pseudoparticles, and the effective pseudoparticle lattices remain invariant under such a transformation. 

The pseudoparticle - pseudofermion unitary transformation is generated by the operator (|81|l . which only shifts the 
pseudoparticle bare-momentum values. Thus, given a pseudoparticle of bare-momentum q, the transformation law 
g — > g of its bare momentum q also defines the transformation law of such a quantum object under the pseudoparticle 
- pseudofermion unitary transformation. For instance, if a specific value of the bare momentum q remains invariant 
under that transformation, then the pseudoparticle carrying bare momentum q is the same quantum object as the 
corresponding pseudofermion. 



A. TRANSFORMATION LAWS OF PSEUDOPARTICLES AT THE BARE-MOMENTUM LIMITING 

VALUES 

The transformation laws of the limiting pseudoparticle bare-momentum values ±qa,,y under the pseudoparticle - 
pseudofermion transformation defined in Eq. (|85|l provide interesting information about the physics described by the 
aiy pseudofermions. For the excited states A and B, the limiting pseudoparticle bare-momentum values /.Qqi/, 
where i = ±1, can be written as. 



L qa,, = L qa„ + I- ^Qaiy ; a = c, s ; i/ = l,2,...; t = ±1 . (124) 

Here q'^^ is the ground-state limiting bare-momentum value given in Eqs. (|39|l and (|40|) . 

We start by considering the case of > pseudoparticle branches. In that case the deviations Aqai, have the 
following form, 



A ^ 



AM„- ^fjy + i/'- l^-i/'l'jAA^a^- ; a = c, s; v = l,2,...] 6 = ±1 , (125) 



where ANai, (and ANau' ) is the deviation in the number of av pseudoparticles (and ai'' pseudoparticles) and AMa 
stands for the deviation in the holon number Mc = [Na — Nco] or spinon number Ms ~ Ncq. These expressions are 
easily obtained by combination of Eqs. (|A3() . HA6|I . (|A7() . and (|A11() of Appendix A. 

In Appendix C we find that at q ^ ^qau and for > the functional defined in Eq. (|73|l equals the following value 



for all excited states A and B, 



X! X! / dq'^a^,a,^>{Lqa^,q')ANa'^'{q') 



= —iAqau; a — c, s ; 1^ — 1,2,...; l — ±1 , (126) 

where the quantity Aqau is given in Eq. H125() . Here it is assumed that the bare momentum q' associated with 
the pseudoparticle bare-momentum distribution function deviations ANa',y' {q') on the right-hand side of Eq. H126|l 
belongs to the domain q' S {—qa'v', +'Za'i/') and can be such that q' — s- ±5^/^' but q' 7^ qa'v'- Interestingly, the value 
of the momentum functional Qau{<l)lL at q = l 

Qaiy is such that it precisely cancels the term l Aqo^jy appearing in the 
limiting bare- momentum expression on the right-hand side of Eq. H124() . Thus at the limiting bare- momentum values 
q — -iiqav given in Eq. I|124l) the value of the momentum functional q{q) defined in Eq. 172|l is independent of the 
value of the deviation Aq^u on the right-hand side of Eq. (|124|l which defines the final excited state. By use of this 
result in Eq. (|72(l we find that. 



Lqau = l^q^av + I'Aqav iqau = i^qau'^ a = c, s ; i/=l,2, t = ±l. (127) 

Thus, we conclude that for v > pseudoparticle branches the limiting momentum values of the corresponding av 
pseudofermion bands equal for all excited states A and B the ground-state limiting bare- momentum values given in Eqs. 
(|39|l and (|40|l . i.e. ztqai, = ^lau- While the width of the bare- momentum Brillouin zone of the av pseudoparticles, 
2qai^, has according to Eqs. H124|) and (|125|l an exotic dependence on the state occupancy configurations, the value 
of the corresponding momentum width of the av pseudofermion Brillouin zone is constant and equals 2qai, = '^Qa^, 
where is the ground-state limiting bare momentum given in Eqs. H39|) and (|40|l . 
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This confirms that the width of the au pseudofermion momentum domain remains unchanged, as in the case of 
the momentum bands of non-interacting particles. This result is consistent with the expressions given in Eq. H77() 
for the rapidity and rapidity-momentum functions, respectively, of the excited states. When expressed in terms of 
the pseudofermion momentum q such functions have the same form both for the ground state and excited states. 
This introduces the requirement that the width of the momentum domain where these functions are defined must be 
the same for the ground state and the excited states. Thus we have just confirmed that this requirement is fulfilled 
through the cancelling associated with the value of the momentum functional H126|l . 

On the other hand, the bare-momentum width of the cO pseudoparticles is independent of the pseudoparticle 
occupancy configurations and reads 2q^Q — 27r/a. In Appendix C we find that the transformation (|85() leaves this 
width invariant and shifts the bare-momentum values q = and q — iq'^Q by the same amount as follows. 



0-^; ±9°o-±'Z°„ + ^, (128) 

where Qcq{0)/L is the q — value of the momentum functional QcQ{q)/L given in Eq. H73(l . Thus, for the cO 
band the shift is the same for these three bare-momentum values. As a result the width of the corresponding 
pseudofermion momentum domain remains unchanged and is given by 2gJ?Q = 27r/a. Note that the value of such 
a shift is a functional which depends on the actual occupancy configurations of the excited states pseudoparticle 
bare- momentum distribution function deviations. 

In general, a ai' pseudoparticle is different from the corresponding av pseudofermion. Below we find that for av 
pseudoparticle branches such that av ^ cO and av ^ si, the only exception is for bare-momentum values q, such that 
g ^ in the particular case when these pseudoparticles are involved in the bare-momentum distribution function 
deviations of zero-momentum excited states. According to the above results concerning the i/ > branches, as 
q -iiqav the excited-state av pseudoparticle limiting bare-momentum values H124() map onto the same momentum 
values q — igj^^, independently of the value of the small deviations ztAg^i, on the right-hand side of Eq. I|124l) . Like the 
corresponding av pseudoparticles 30], as g — > ig^^ = igj^^^ the av pseudofermion and the rotated av pseudofermion 
become the same quantum object. In that limit, the energy tcv{q) ^ 2fj,v + e°j,(g) (and £31/(9) = + £^1/(9)) of 

the composite cv pseudofermion (and sv pseudofermion) is additive in the —1/2 Yang holon energy 2fj, (and —1/2 HL 
spinon energy 2^oH) because e^^{q) (and £%{q) ^ 0) as q ^ ±qw — ±9ci/ (and q ±qsi^ = i^D- It follows 
that in that limit the v —1/2 holons (a = c) or —1/2 spinous {a = s) contained in the composite av pseudofermion 



w ps 

acquire the same localized character as the —1/2 Yang holons (a = c) or —1/2 HL spinous {a = s) |33, respectively. 



B. TRANSFORMATION LAWS OF FERMI-POINT PSEUDOPARTICLES 



Let us next consider excitations corresponding to J-CPHS ground states whose Fermi points are defined in Eqs. 
(|69|l - l|71|l . Our goal is to find out how the Fermi bare momentum qFcu,L given in these equations transforms under the 
pseudoparticle - pseudofermion transformation 185(1 . In this case, such a transformation can formally be written as, 

where according to Eq. (|71|l . Aqpa^.i — t [2'^/L] lS.Nav,L- In Appendix C, we find that the deviation lS.qpa,j,u on the 
right-hand side of Eq. (|129|l is given by, 

Ag>„., = . ^ AAT/^^ , ; AN^, . = + ^ ■ (130) 

Here 



CL.a'.'AJcV, (131) 
and 



A7V„, 



^jL - E E 



a.'—c,s u'—5^ 



CL,a'.' ='^o,o'5.,.'+ E (0^'C'.'(4a.,^''zLv'); J=0, 1 

= ±1 



(132) 
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where two-pseudofermion phase shift given in Eq. ifSHl . 

The pseudoparticle - pseudofermion transformation (|85|l . refers to a subspace spanned by states of finite energy uj. 
Let us show that in the specific hmit of small energy w — > that unitary transformation is directly related to the 
primary-field conformal dimensions of the ID Hubbard model two-component conformal field theory |^|23,|23|- We 
consider the particular case of excited states with finite occupancies in the cO and si pseudoparticle bands only. For 
electronic densities and spin densities such that < n < 1/a and < m < n, respectively, this corresponds to the 
low-energy Hilbert subspace considered in the studies of Refs. [23.l23j . For simplicity, in this specific case we use the 
notation of these references such that c = c, and s = s, 1. The indices av of all quantities are replaced by the index 
a such that a = c, s. In the present low-energy Hilbert subspace the above general expressions simplify to. 



where 



and 



(134) 



E 

a'—c,s 



.0 
Sa, 



AJl 



(135) 



Here 



^la'=^".«'+ E (136) 
i,'=±l 

Let us combine Eqs. H134|l and H135() and rewrite the number AA'^ ^ as. 



[E C 



+ L 



E ^ia'AJ. 



(137) 



Based on the form H136I) of the parameters S^l^ ^, one can by manipulation of the integral equations given in Appendix 
B, define these parameters in terms of related integral equations. Importantly, comparison of these latter equations 
and associated expressions ()136|l and (|137|l with the results of Refs. ,2^i2^J reveals that the number AA^^ ^ given in 
Eq. (|137|l is such that, 



[AAT/ = 2 AJ, = [ 5: „^ ^ + ^ E A J.,] ' , (138) 



where the quantity A^ is the a conformal dimension of the primary fields |20l |21| associated with the a = c and 
a = s excitation branches. Moreover, the parameters ^/ and ^JJ ^, are entries of the transpose of the dressed charge 
matrix and of the inverse of the transpose of the dressed charge matrix, respectively, of the low-energy two-component 
ID Hubbard model conformal field theory [2ll2ll23|. 

These results confirm that for the low-energy subspace defined above and for electronic densities and spin densities 
such that < n < 1/a and < m < n, respectively, the conformal invariance of the ID Hubbard model is directly 
related to the pseudoparticle - pseudofermion unitary transformation (|85|l . For these values of the electronic and spin 
densities the low-energy excited states are described by pseudoparticle bare-momentum occupancy configurations 
such that there are both pseudoparticle and pseudoparticle holes in the c = c, and s = s, 1 pseudoparticle bands. 
However, the pseudoparticle - pseudofermion unitary transformation is more general and refers to all values of energy 
associated with few-electron excitations. 

The above results provide interesting information about the transformation laws of cO pseudoparticles (and si 
pseudoparticles) of bare momentum values at the Fermi points, q — i2kF + i- ^ AA^cO, t (and q — l kpi + i- ^ ATV^i^ J 
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where l = ±1. Under the pseudoparticle - pseudofermion transformation (|85|l . these pseudoparticles are mapped onto 
cO pseudofermions (and si pseudofermions) with momcntmn values at the Fermi points q — L2kp + l Aj; (and 

q — ikpi + L ^^72 A^) where A^ (and A!,) is the conformal dimension of the cO (and s) excitation branch primary 
field [23, 122, 123 ■ We then conclude that these conformal dimensions are such that the quantities A^ are nothing 
but deviations in the values of cO (a — c) and si (a = s) pseudofermion momentum Fermi points resulting from low- 
energy excitations. (We note that the positive and negative value of the root yj2 A^ refers to pseudofermion creation 
and annihilation, respectively.) Conformal-field theory can be used to evaluate expressions for low-energy few-electron 
spectral functions |2ll . |2^ . We recall that our general pseudofermion description provides the momentum deviation 
values resulting from finite-energy excitations for all values of pseudofermion momentum and for all pseudofermion 
branches. Then it is not unexpected that the general pseudofermion description introduced in this paper can be used 
in the evaluation of few-electron spectral function expressions for all values of energy L> lQ.-> ^ ■ 



C. INVARIANCE UNDER THE PSEUDOPARTICLE - PSEUDOFERMION UNITARY 

TRANSFORMATION 

According to the results of Appendix A, there are no ciy pseudoparticles and sv pseudoparticles belonging to 1/ > 
and v > 1 branches, respectively, in the ground state. This is consistent with the Fermi values of Eq. I|()8|) . such that 
iFau ~ ^ ^'^^ these av pseudoparticle branches. In Appendix C, it is shown that the momentum functional Qai,{Q)/ L 
given in Eq. (|75|l vanishes at q^-,^^ = for these av pseudoparticle branches, 

9l^=Y, p'"' dq<i>^,,^,,,{0,q')AN^,,,{q)^Q, (139) 

a'—c,s u' — 1 — 5^1 ^ ~^a'v' 

provided that the excited states are partial J-CPHS ground states and the bare-momentum distribution function 
deviations on the right-hand side of Eq. H139|) are for the cO and si pseudoparticle branches such that, 

AAr,o(g) = AA^co(-g) ; AiV,i(q) = AiV,i(-g). (140) 

Such excited states have vanishing momentum. We recall that partial J-CPHS ground states are states that have 
J-CPHS ground-state occupancy configurations for all the av pseudoparticle branches other than the cO and si 
pseudoparticle branches. It follows that for the former branches the vanishing pseudoparticle Fermi points, 



qPav, L = t q^au + ^qPau, ' ^ ^ + ^N^u, t , (141) 

remain invariant under the pseudoparticle - pseudofermion transformation H85|l . or 



qPay^L ^ qPau^L — qPav.L ■ (142) 

This means that a au pseudoparticle belonging to a branch other than the cO and si branches, carrying bare 
momentum q — s- 0, and involved in excitation processes associated with zero- momentum excited states whose cO and 
si bare- momentum distribution function deviations obey Eq. (I140|l remains invariant under the pseudoparticle - 
pseudofermion transformation 1)8 5() . It follows that in this limiting case the q ^ av pseudoparticle is the same 
quantum object as the corresponding q — s- av pseudofermion and then, 



q q = q; as q ^ . (143) 

This invariance occurring as g — > for av pseudoparticles belonging to branches such that av ^ cO and av ^ si, 
is associated with the free, non-interacting, and delocalized character of these quantum objects. Indeed, in this limit 
the residual interactions vanish and these pseudoparticles are free, non interacting, and delocalized quantum objects 
provided that the excitation processes they are involved in correspond to excited states of zero momentum. On the 
other hand, according to the discussions presented in Sec. HI, as q approaches the limiting bare-momentum values 
()124)l . q ±qa^, the av pseudoparticles belonging to branches other than the cO pseudoparticle branch become 
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non-interacting and localized. Such & q ^ iqav non-interacting and localized behavior is associated with another 
invariance: As g — > iqau these av pseudoparticles become invariant under the electron - rotatcd-electron unitary 
transformation |3nl l8l|. Thus, as the bare momentum q approaches the limiting values iqau the pseudoparticle 
and the rotated av pseudoparticle become the same quantum object. In contrast, according to Eq. (fHTjl in that 
limit, these av pseudoparticles are not invariant under the pseudoparticle - pseudofermion unitary transformation. 
Their non-interacting and localized behavior is achieved by a cancellation of the term t Aqai, appearing in the limiting 
bare-momentum expression on the right-hand side of Eq. H124|l by the momentum functional Qau^i^qau) I L of Eq. 
(|126|l . Moreover, for other bare momentum values q such that g ^ and q ^ ±qau the av pseudoparticle is different 
both from the av pseudofermion and rotated av pseudoparticle. 

We thus conclude that as g — > a av pseudoparticle involved in deviations associated with zero-momentum excited 
states becomes free, non-interacting, delocalized, and the same quantum object as the corresponding av pseud- 
ofermion. As the limiting values are approached and thus q — > iqo,^, the av pseudoparticle becomes non- interacting, 
localized, and the same quantum object as the corresponding rotated av pseudoparticle. In this latter case, the energy 
of this composite object becomes additive in the energy of its v —1/2 holons (a — c) or v —1/2 spinons {a — s). 
Note that one of the effects of increasing the bare- momentum absolute value from \q\ = to \q\ = qav is to enhance 
the localization character of the av pseudoparticles. If a pseudoparticle has bare-momentum q and q + Aq for an 
initial ground state and an excited state, respectively, then the limiting free, non-interacting, and delocalized character 
and non-interacting and localized character are achieved when Qav{q)/L = and Qav{q)/L = — Ag, respectively, 
where Qau{q)lL is the momentum functional (|73|) . For other values of that momentum functional the pseudoparticle 
is interacting. The pseudoparticles are never invariant under both the pseudoparticle - pseudofermion and pseu- 
doparticle - rotated-pseudoparticle unitary transformations. This is because invariance under the pseudoparticle - 
pseudofermion unitary transformation is associated with a free, non-interacting, and delocalized character, whereas 
invariance under the pseudoparticle - rotated-pseudoparticle unitary transformation corresponds to a non-interacting 
and localized character. Therefore the impossibility of simultaneous invariance under these two transformations is 
related to Heisenberg's uncertainty relation. 

Another interesting example of invariance under the pseudoparticle - pseudofermion transformation (|85|l concerns 
the bare-momentum difference {q — q') when q and q' differ by a finite multiple of 2-k/L. This case is of physical 
importance because (g — q') corresponds to the momentum of a pseudoparticle - pseudoparticle hole excitation in the 
cO or si band. This type of process is associated with the excited states B and plays a key role in the few-electron 
spectral- weight distribution |34j . In this case, we find in Appendix C that in the thermodynamic limit, {q — q') is 
indeed invariant under the pseudoparticle - pseudofermion transformation and thus, 

{q'q')~^{q-q')^{q-q')- (144) 

VIII. CONCLUDING REMARKS 

In this paper we introduced a pseudofermion operational description for the ID Hubbard model. We found that 
in the thermodynamic limit the wave function of excited states belonging to the Hilbert subspace of the ground- 
state normal-ordered ID Hubbard model associated with few-electron excitations factorizes for all values of U /t. 
This factorization results from the non-interacting character of the pseudofermions whose oc cup ancy configurations 
describe these excited states. While the pseudoparticle description studied in Refs. [33, |^ corresponds to the 
whole Hilbert space, the pseudofermions are defined in the Hilbert subspace associated with few-electron excitations. 
In such a subspace, the pseudofermions are related to the pseudoparticles by a unitary transformation. We classified 
the statistics of the latter quantum objects according to a generalized Pauli principle [36j |. 

Our study included the introduction of the pseudoparticle - pseudofermion unitary transformation and of an operator 
algebra for both the pseudoparticles and pseudofermions. Such a transformation removes the residual interactions 
of the av pseudoparticles, which are mapped onto the non-interacting av pseudofermions. The av pseudoparticle 
residual interactions are cancelled by the momentum functional Qa^v{q) / L of Eq. H73I) . The information recorded in 
the pseudoparticle interactions is transferred over to the momentum two-pseudofcrmion phase shifts of that functional. 
These phase shifts control the few-electron spectral properties through the same unconventional momentum functional. 

We introduced creation and annihilation operators for both the pseudoparticles and pseudofermions and derived the 
anticommutation relations of these pseudofermion operators. We also studied and discussed the transformation laws 
of the pseudoparticles under the pseudoparticle - pseudofermion transformation. This study included the discussion 
of the physics behind both these transformation laws and the invariance under the pseudoparticle - pseudofermion 
and pseudoparticle - rotated-pseudoparticle unitary transformations of the pseudoparticles for some specific bare- 
momentum values. Invariance under the former (and the latter) transformation is associated with a free, non- 
interacting, and delocalized character (and non-interacting and localized character) for these quantum objects. Thus 
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as a result of Heisenberg's uncertainty relation, these two invariances never occur simultaneously. We also find that 
under the pseudoparticle - pseudofermion unitary transformation, the cO and si pseudoparticles of bare momentum 
at the Fermi points are mapped onto corresponding pseudofermions whose momentum Fermi points expressions are 
related to the conformal dimensions of the two-component conformal- field theory primary fields |23,|23|- It follows that 
in the limit of low energy the conformal invariance of the ID Hubbard model is related to the general pseudoparticle 
- pseudofermion unitary transformation. 

The pseudofermion algebra introduced in this paper is used elsewhere in the evaluation of finite-energy few-electron 
spectral- function expressions 0, |3 ■ Fortunately, as a consequence of the above wave function factorization, the 
few-electron spectral functions can be written as a convolution of pseudofermion, —1/2 Yang holon, and/or —1/2 HL 
spinon spectral functions for all values of energy and on-site repulsion U [s^l • 



Acknowledgments 

I thank Karlo Penc for many useful and illuminating discussions concerning the issues studied in this paper. I 
also thank Daniel Bozi, Ralph Claessen, Francisco (Paco) Guinea, Katrina Hibberd, Luis Miguel Martelo, Pedro 
Sacramento, and Joao Lopes dos Santos for stimulating discussions. 



APPENDIX A: PSEUDOPARTICLE REPRESENTATION, EFFECTIVE PSEUDOPARTICLE LATTICE, 

AND GROUND-STATE QUANTITIES 

In this Appendix, we summarize some aspects of the pseudoparticle, description which are needed for the studies of 
this paper. This includes a brief summary of the basic pseudoparticle properties, introduction to the pseudoparticle 
bare momentum jsj, local pseudoparticle concept, and effective lattices [sEj- In addition, we provide the ground-state 
pseudoparticle bare-momentum distribution functions, as well as other ground-state quantities. 

According to the results of Ref. I.Sfl . the at^ pseudoparticle bare-momentum description obtained naturally from the 
Bethe-ansatz solution in Refs. |22ll23. 29, 31] is related by Fourier transform to a local av pseudoparticle description 
in terms of spatial coordinates of an effective av lattice. These concepts are needed and useful for both the operational 
pseudoparticle and pseudofermion representations introduced in this paper. 

The aiy pseudoparticles carry bare momentum q. This is the continuum bare momentum associated with the 
discrete bare-momentum values qj such that. 

These discrete values read [Slf . 

9.-^^r' (A2) 

where the /""^ numbers are integers or half-odd integers |30l IsH . The index number j can have the values j = 
1,2, ...,N*^. Here the number N*^ equals the number of discrete bare-momentum values in the aiy pseudoparticle 
band and is such that, 

- N^, + Nt , (A3) 

where Nai, is the number of aiy pseudoparticles and iV^^ is the number of aiy' pseudoparticle holes. The latter number 
is given by, 

oo oo 

iV^'o ^Na-N,o; N!i, = 2S^ + 2 ^ {ly' - v)N^,, =L^ + 2 ^ {v' - v)N^,, ; a = c, s , > , (A4) 

where Sc and Ss are the values of ?7-spin and spin, respectively, which read 

^ oo ^ oo 
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Combination of Eqs. I|A4|I and ljA5|l leads to, 

For ai' 7^ cO (and ai' = cO) the numbers Hj"^ on the right-hand side of Eq. ljA2|l are integers (half-odd integers), 
if N*^, (and Na/2 — J2a=c s^'uL^ ^au) is odd (even). This non-perturbative shake-up effect is related to the so 
called orthogonality catastrophe |[41| . It reveals that the values of the available bare- momentum values qj of each av 
pseudoparticle band might be different for different CPHS ensemble subspaces. 

The holons, spinous, ar id p seudoparticles are related to the electrons through the rotated electrons, as discussed in 
Sec. Ill and in Ref. |^ • The unitary transformation which maps electrons onto rotated electrons is defined in 
Sec. II. Since the momentum operator is invariant under such a transformation, also the lattice occupied by rotated 
electrons has the same lattice constant a and length L as the original electronic lattice. In Ref. 35] it was shown 
that there is an effective av pseudoparticle lattice for each av pseudoparticle branch. 

In contrast to the case of a non-interacting band, the value of the number of available discrete bare-momentum 
values N^jj defined in Eqs. (IA3|I . I)A4|I . and (|A6|) might be different for different CPHS ensemble subspaces. This 
value was identified in Ref. 35] with the number of sites of the corresponding effective av pseudoparticle lattice. This 
lattice has the same length L — N*^ a^i, as the original real-space lattice, where 

= = (A7) 



is the effective av lattice constant. Since A^*^ is different for different CPHS ensemble subspaces, it follows from Eq. 
(jA7|l that the value of the corresponding effective av lattice constant a^v also changes. The spatial coordinates of 
the effective av pseudoparticle lattice are Xj = aavj where j = 1, 2, 3, A^*^. The value of A^*^ defined in Eqs. 
(|A3(I - (|A6(I determines the corresponding limiting bare-momentum values of the olv pseudoparticle Brillouin zone. 
In the case of the cO pseudoparticle band, the discrete bare-momentum values qj belong to the following range, 

< qj < qto , (A8) 

Here, 

.cV-.c-„-^[l-^], (A9) 

for Nco even and 

C.^-^[l-^]. (AlO) 

for NcO odd. Note that since N*q = Na, the lattice constant of the effective cO lattice equals the lattice constant a of 
the original electronic lattice, that is Oco = a. Thus, for the cO band the lattice constant Oco has the same value for 
all CPHS ensemble subspaces. 

For the av pseudoparticle bands such that v > 0, the discrete bare- momentum values qj are distributed symmet- 
rically relative to zero, and are such that \qj\ < q^i,. The two band momenta ±qai, are the limiting bare-momentum 
values associated with the limits of the av pseudoparticle Brillouin zone. The limiting bare- momentum Qai, reads, 

= (AH) 

where the lattice constant Oav > a is given in Eq. (|A7ll . 

We close this Appendix by providing some useful ground-state quantities. Following the results of Refs. l30l]. 
the ground-state pseudoparticle bare-momentum distribution functions and —1/2 Yang holon and —1/2 HL spinon 
numbers read, 

A^co(9) = e(g+,o-(?), 0<(z<g+; iV°o («) = © ^ ^fco) > 9co < 9 < ; (A12) 
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N°,iq) = e(^qFsl-q) , < q < Qsl ] N°,{q) = e(^q + qpsl) , -qsl<q<0; (A13) 

^Liq) = -1/2^0; ~qa^<q<qau\ a = C,S; iy>l+Sa,s- (A14) 

Except for l/L corrections, the Fermi values on the right-hand side of Eq. (jA12|l are given by q$^Q = ±qp^Q, where 
the ground-state values qp^^ and g^^g are given in Eq. . Here we have ignored terms of order of 1 /L and used the 
expressions kpa- ~ ttUo- and 2kp = nn. The expressions of the Fermi values including l/L contributions are provided 
in Appendix C of Ref. ,31j . 

For the ground-state CPHS ensemble subspace the expression of the number N*^^ given in Eqs. (|A3(I . (|A4(I . and 
(jA6|) simplifies. Let us denote by N^'* the ground-state value of the numbers N*^. These numbers read, 

^cV = (^a-iV); N^i*^Nt, iV°'* = (iVt - iVi) , i.>l, (A15) 

whereas N^!^* = N*^ is given by 7V*q = Na for the whole Hilbcrt space. 

In addition to the —1/2 Yang holon and —1/2 HL spinon numbers already given in Eq. (|A14(I . the ground state 
belongs to a CPHS ensemble subspace with the following values for the pseudoparticle, ±1/2 holon, and ±1/2 spinon 
numbers, 

M0_i/2-0; A4%i/2 = L^_+i/2 = 7V,-7V; iV," = TV ; = Q , (A16) 
in the charge sector and 

A4Vi/2 = = ; A4"+i/2 = a^t; = o, v>2- lI+,/^ = - n^, (A17) 

in the spin sector. 

Finally, note that for > the number ljA3|l can be expressed in terms of the ground-state values ()A15|I as follows, 

oo 

Ku ^ + ALo, + 2 J2 W-v)No..'\ a = c,s, i/ > . (A18) 
Moreover, for aiy — si the pseudofermion hole number defined in Eqs. (|A4() and HA6|I can be written as, 

oo 

^N,o~2j2 Nsw ■ (A19) 

iy' = l 

APPENDIX B: THE TWO-PSEUDOFERMION PHASE SHIFTS (r,r') 

Here we provide the set of integral equations which define the two-pseudofermion phase shifts ^af, a'v' [f, f') on the 
right-hand side of Eqs. lEOIl-linB, GH), and (jSljl . 

Let us start by introducing the following Fermi surface parameters, 

n 4i sin Q At B , , 

where the parameters Q and B are defined in Eq. (|45|l . 

In order to derive the integral equations which define the two-pseudofermion phase shifts ^av, a'v> {f, r'), we first use 
in Eqs. (|19l) - (|21|l the bare-momentum distribution functions of the general form (|59|) as well as the rapidity functional 
expressions given in Eq. (|77|l and rapidity-momentum functional expression provided in Eq. |(77J). Such functionals 
are written in terms of corresponding ground-state rapidity and rapidity-momentum functions whose argument is the 
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momentum functional. The expression of this functional is given in Eq. H72I) with the momentum Qaviq)/ L provided 
in Eq. (|73|l . Expansion of the obtained equations up to first order in the bare-momentum distribution function 
deviations on the right-hand side of Eqs. (|5H|l and l(7!!jl leads to expression (fZ^ with the two-pseudofermion phase 
shift q';/' (r, r') uniquely defined by the integral equations given below. A first group of two-pseudofermion phase 
shifts obey integral equations by their own. These equations read, 



$si.co(r,r') ---arctan(r-r')+ / ' dr" G(r, r") ^ai, cO (/', r') , (B2) 



TT 



1 



arctan ( — 



^si,cu{r,r')^-— dr" / / + / dr" G(r, r") <i>«i, (r", r') , (B3) 



and 



* / ,x 5i V fr-r'\ (l-6iy)( /r-r' 
9si,su{r,r) = arctan — - — H ^ — <^ arctan 



arctan ( ) \ 



1 



dr"- 



arctan 



+ dr"G{r,r")^,,.si{r",r') . 



(B4) 



Here the kernel G{r, r') is given by, 



G(r,r') 



2tt 



l + ((r-r')/2)2 



(B5) 



where 



t(r) = -[arctan(. + r«)-arctan(r-r°)] , 



(B6) 



and 



lir) = - [HI + {r + rlf ) ln(l + (r - r^f)] 



(B7) 



The kernel defined in Eqs. ljB5|) - (|B7|) was first introduced in Ref. '23] within the low-energy two-component c = cO 
and s = si pseudoparticle theory studied in that reference. 

A second group of two-pseudofermion phase shifts are expressed in terms of the basic functions given in Eqs. 
dHa-dH as follows, 



$cO,cO {r, r') 



dr' 



l + {r- r" f 



(B8) 



- 1 { r — r' 
^cQ,cu {r, r') = arctan 

TT \ V 



dr' 



H- (r - r"f ' 



(B9) 



and 



^co,.. M = -iarctanf^) -f i T d/'^^i^fi^ . (BIO) 



TT \ V / TT 
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Finally, the remaining two-pseudofermion phase shifts can be expressed either in terms of the functions ljB8l) - (|B10|) 
only, 

^0 O - -arctan(-- j - - dr ^ ^r_^y^ , (BH) 
$c.,c.' (r,/) = f e.,.,(r - r') - - T ^^''^T^^ , (B12) 

and 

or both in terms of the basic functions ljB2() - (|B4|l and of the phase shifts (|B8|l - (|B10(l . 

TT J_rO + J-rO 27r 

and 

In the above two-pseudofermion phase shift expressions the functions and 0|)'j^, (a;) read, 

0!y,!y'(a;) = ^iy,;,'|2arctan(^) + y^ 4arctan(|y) | + (1 - ^/)|2 arctanf ^ ^ ^ ) 



2 arctan 



and 



'^''^'^''^ " -''^''-'Un +(^)2i +2^zh +(|.)2i/ + ^^ '^-.-'^1 



2 
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i^ + ^')[^ + ii;fp?] t[ (l^-^'l + 20[l + (^;^T^)^] 

respectively. Note that the latter function is the x derivative of the function defined in Eq. ljB17|l . 

In spite of the different notation and except for simplifications introduced here as a result of some integrations 
performed analytically, the integral equations H5U|I - H51() are equivalent to the system of coupled integral equations 
(B30)-(B40) of Ref. 
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The bare-momentum two-pseudofermion phase shifts are defined by Eq. H74I) . in terms of the above phase shifts 
^aiy,a'v'{r, r') associated with the integral equations l|B2l) - (|B15|) . In apphcations of the pseudofermion description 
to the study of few-electron spectral properties the phase shifts $si,si(9, '(),J^s1.m(q q'), $cO,co(Q) <i')i '^cO,si{q, q'), 
^co, ci{q, q'), and ^si, ci{q, q') with q at the Fermi points play a major role By manipulation of the above integral 
equations, we find that in the limits of zero spin density m — > and U /t 0, these two-pseudofermion phase shifts 
with q at the Fermi points and the second bare-momentum denoted by q are given by, 



^si.siiikp, q) 



2^2 



qj^Lkp; 



2V2 



q = bkp , L — ±1 . 



(B19) 



^si,ca{i'kF.q) = ^{-ie(2/cj.-|g|)-f sgn(g)e(|g|-2fcF)|, q^i2kF\ 

L 

q = L 2kF , i = ±1 , 



2V2 ' 



(B20) 



,,o(t2fcj.,g) = ^^^~ie{2kF-\q\)+Sgn{q)Q{\q\~2kF)} , q^i2ki 

r 3 1 

- ij , q^ L2kF , t = ±1 , 



(B21) 



^ CO., slit '2k F, q) 



2V2 
sgn(g) 
2V2 



q^ Lkp; 
, q ~ ikF , i = ±1 



(B22) 



and 



*cO,ci(t2fcF, q) 



sgn(g) 
V2 



q^O; 



^si.ciii-kF, q) = 0, i = ±1 , 



(B23) 



(B24) 



respectively, where 6(a;) = 1 for x > 0, and 6(x) = for a; < 0. 

For m ^ and U/t 00 the phase-shift expression ljB19|l for ^sLsiii-kp, q) remains valid, whereas the two- 
pseudofermion phase shifts $si,co(g, q'), ^co.coiq, q'), ^•si,ci(g, q'), '^cO.,si{q, q'), and ^co.cite, q') with q at the 
Fermi points and the second bare-momentum denoted by q read, 



^si,ca{i^kF, g) = - 



2V2 



t = ±1. 



(B25) 



$cO,co('-2fc_F, q) = ^si^ciiikF, q) =0, t = ±1 , 



(B26) 



$cO,si(t2/cF, q) = -rr- , 
Akp 



\q\ <kF , t = ±1 , 



(B27) 



and 



*cO,ci(t2fcF, q) 



2[n- 2kF] 



\q\<7r- 2kF , 



±1. 



(B28) 
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APPENDIX C: TRANSFORMATION LAWS UNDER THE PSEUDOPARTICLE - PSEUDOFERMION 

SUBSPACE UNITARY ROTATION 



Here we study the momentum functional Qai'{(i)/L given in Eq. I|73(l for specific values of q. In some cases we 
consider general excited states obeying relations H7()|) . whereas in other cases we consider particular cases of such 
general excited states. 

We start by confirming that the second term on the right-hand side of Eq. H79(l is of order. Let us consider 

a more general situation and also confirm the validity of Eq. H144|l and show that {q — q') is invariant under the 
pseudoparticle - pseudofermion transformation when q = q' + Nph^i: / L] and Nph — ±1, ±2,... is a finite integer 
number. Equation (|79|l corresponds to the particular case when Np^ = 1. Thus we want to show that the quantity 
[Qau{<l + ^ph^) — Qau{Q)]/L involving the functional (f7^ is of {l/L)^ order. By expressing this quantity in terms 
of the derivative dQa^{q)/dq one finds, 



Jl^r>.^. (CI) 

Analysis of the form of the derivative dQau{q) / dq reveals that it is of order. This is the confirmation that the 

momentum contribution IjCip is of second order in 1/L. Thus, since within the present large-L pseudofermion descrip- 
tion bare-momentum and momentum contributions vanish, {q — q') is indeed invariant under the pseudoparticle 
- pseudofermion transformation. 

The remaining of this Appendix is complementary to the studies of Sec. VI about the transformation laws of the 
pseudoparticle - pseudofermion subspace unitary rotation. Next we consider that q = ±qau for v > Q pseudoparticle 
branches. Our goal is the study of the momentum functional. 



J2 ^ dq' ^au.,a'y'{iqay,q')^Na'u'{q')] a = c,s; v>0; t = ±l, (C2) 



Q'— c,s u' —1 — 6 



and confirm that it is given by expression H126|l where the deviation Agci^ is provided in Eq. 1)125(1 . It is assumed that 
the bare momentum q' belongs to the domain q' e (— gav, +qa'u') and can be such that q' — > ±qa'u' but q' ^ qa'v' ■ 
First we note that according to Eq. H42I) . the two-pseudofermion phase shift ^a_v.a'v\^qavT<i^ expressed in terms 
of the pseudoparticle bare momentum on the right-hand side of Eq. (|C2|I is such that, 



$a.. (±ga., = aV ( ±00, ) , (C3) 



where the two-pseudofermion phase shift ^av.a'v' is defined in Appendix B. In order to achieve this result we 

used that A^^(±g„i,) = ±oo, as given in Eq. W^ . 

By manipulation of the integral equations given in Appendix B, we find that all the two-pseudofermion phase shifts 
of form l(C3|l vanish except the following ones, 

^at/.co(±gai.,g') ^± ^^"''' ^ ^"'"^ ; $^^^^^,(±q„^,g') = zp^^l^l ± ^; a = c, s ; v,v'>^. (C4) 

Use of Eq. HC4|I on the right-hand side of Eq. HC2|I leads to. 



L = ^2 



OO 

ANc,^-[Sa,c-Sa,s]ANco-^(iy + i^' -\i^-i^'\)ANc,.']^T^q»i^; a = c,s; i/ > , (C5) 



where we also used expression (I125|l . Finally, note that expression (|C5p is equivalent to Eq. (|126|l . 
Let us now consider the quantity, 



L 



E 



c.s iy' — l — 5^ 



dq' $cO,a'.'(gc1).'7') AiVaV'(g'), 



(C6) 
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and confirm that it obeys relation (|128|l . This is easily confirmed by noting that according to Eq. (|42|l the two- 



pseudofermion phase shifts ^cO,a'i^'iqfo, l') on the right-hand side of Eq. HC2|I are such that, 



0, 



A°o('7'; 



(C7) 



where the two-pseudofermion phase shifts ^co.a'u' t"') are defined in Appendix B. The result (|C7p implies the 
validity of Eq. (fngjl . 

Our next task involves evaluation of the momentum functional Qaiy{Q)/L given in Eq. I|73|) for q = Lqpai, in 
the particular case when the excited states associated with the deviations ANa',^'{q') on the right-hand side of that 
equation are J-CPHS ground states. The Fermi bare momentum values that limit the compact pseudoparticle bare- 
momentum occupancy configurations of these states are given in Eqs. ifTUI) and (|7T|l . While the above expressions 
derived in this Appendix refer to general excited states whose deviations obey relations (|7^ . the following results 
are valid only for the particular case when these excited states are J-CPHS ground-states. Our goal is to arrive to 
expressions (|130f) - (ll32f) . This implies evaluation of the following quantity, 



2tt ^ 



(C8) 



where 



L 



E E 



dq' '^au,a'y'{iqFav:q') ANa'„'{q') , 



(C9) 



and the bare-momentum distribution function deviations ANa'^'{q') correspond to final J-CPHS ground states. 
The phase shifts ^ai^^a'i/'{q,q') and ^aiy.a'u' {f^r') have the following property. 



qV (<7, q') = -^au, qV [~q, ~q') {r, r') = -^ai^. a'y' -r') 



(CIO) 



This symmetry is found by analysis of the integral equations given in Appendix B, which define the two-pseudofermion 
phase shifts $qi/, aV associated with the phase shifts $01^. qV (^i through Eq. (|74(l . Once the ground-state 

rapidity functions that appear on the right-hand side of Eq. (|74|) in the argument of the two-pseudofermion phase 
shifts ^av.a'v' are odd functions of the bare momentum, the second relation given in Eq. (jClOp implies the 

validity of the first relation of the same equation. 

Since the pseudoparticle bare-momentum distribution function deviations of the J-CPHS ground-states include 
creation or annihilation of pseudoparticles in the vicinity of the Fermi points only, we can replace q' by sgn(g') qp^i^i 
in the argument of the phase shift ^av, a' v' {i- q%au^ q') '^'^ right-hand side of Eq. I|('9(l . This leads to. 



27r ■' 

AqFc,,^, = L—AN^,,,+ EE/ 0^g'AiVaV'(g')*a.,aV'('•9Fa^.:Sgn(g')gLv')• (dl) 

a' —C. S U'— 1 — 5^1 ^ ^^a'w' 



By performing the q' integrations and using Eq. IjClOp we arrive to. 



AqFav,i, = —j- ^ ^ ^ yi'" 5a,a' 5,,^l" + l''^av,a'v'{q%'ai'T''''" q^'a'u') 

a'—c, s f' — 1 — 6^1 ^ — 

After performing the l" summation and using again Eq. (|(]10|l we find, 



AJ„v']. (C12) 



a'—c. s u' — 1—5^1 



t' = ±l 



27r 



E/ E/ ^a,a'Su,v'+ t-' '^av.a'u'iq^'auT l-' q^'a'u') 



a — c, s u 



2 

AJf^'i/i . 



i,'=±i 



(C13) 
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Note that since Ag^ai., = L ^AiV4, , this result is equivalent to Eqs. ltnni) - ((n^ . 

Finally, let us confirm the validity of Eq. H139|l and show that the momentum functional Qai, (0) / L given in Eq. 
vanishes for av pseudoparticle branches such that a,v ^ cQ and a,u ^ si provided that the deviations AiVco(<;) and 
AA^si('z) on the right-hand side of that equation obey Eq. H14()(l . In this case, the excited states are partial J-CPHS 
ground states. These are states that have J-CPHS ground-state occupancy configurations for av pseudoparticles such 
that a,v ^ d) and o.^v ^ si. 

There are no av pseudoparticles belonging to branches other than cO and si in the initial ground state. It follows that 
excited states obeying relations ((7^ and having J-CPHS ground-state occupancy configurations for av pseudoparticles 
such that a,v ^ and a^v ^ s\ have only a vanishing density of these quantum objects. Such J-CPHS ground 
state occupancy configurations correspond to compact bare-momentum domains of vanishing width centered at g = 
^Fav ~ "^^^ generation from the ground state of these excited states includes creation of av pseudoparticles such 
that a,v ^ d) and a, 7^ si for bare momentum g — > 0. A simple and useful example, is the creation of a single av 
pseudoparticle at bare momentum q = (fp,^^ = 0. To start with, we assume that the deviations ANco{q) and ANsi{q) 
are small but have arbitrary values, whereas the deviations ANai^{q) of the remaining excited-state occupied av 
pseudoparticle branches correspond to J-CPHS ground state occupancy configurations associated with Fermi points 
of the form given in Eq. (|141|l . 

Use of the integral equations provided in Appendix B leads in the case of such partial J-CPHS ground states to, 

$Qi/, qv(Oi'''0) = 0; V, v' > for a, a' = c; v, v' > 1 for a, a' ^ s . (C14) 
Use of Eq. (ICTl m expression ifT^ leads to, 

CL,qV' = ^a,a' Vi'' ; V, v' >0 for a, a' =c; v, v' > 1 for a, a' = s ; j = 0, 1 . (C15) 
Moreover, it follows from Eq. (jClOp that, 

$a.,co(0,g') = -*a.,co(0,-<z'); <^a., sl{0, q') ^ ^^c.. sl{0, -q') ] y>5^..s. (C16) 

Let us next consider the same type of excited states but with small deviations ANca{q) and ANsi{q) obeying Eq. 
(|14U|I . These excited states have vanishing momentum and their cO and si current number deviations AJco and AJsi, 
respectively, vanish, i.e. AJco = AJji — 0. In this case it follows from Eqs. I|C14|) - (|C15|I that. 



dg'$a.,co(0,g')AiV,o(g')+ dg' si(0, g') AiV,i(<z') 
+ T E E E *-.aV'(0,^'0)^i^+ ^'$o.,a'.'(0,t'0)AJ„v']=0. (C17) 

a'— c, s u' — 1+5^1 ^ z,'— itl t'— ±1 

This result confirms the validity of Eq. H139|) . 
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